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ABSTRACT 


Recent  applications  of  laminated  composites  to  ship  structural 
components  have  highlighted  both  the  analytical  complexities  in¬ 
volved  and  the  need  for  providing  the  structural  designer  and 
analyst  with  appropriate  laminate  analysis  methods.  To  this  end, 
three  types  of  analytical  methods,  traditionally  applied  to 
metallic  ship  structural  components,  have  been  extended  to  the 
analysis  of  laminated  composites;  procedures  for  their  implemen¬ 
tation  are  discussed.  These  methods  include  (a)  strength-of- 
materials  techniques  for  laminated  beam  and  column  members,  (b) 
classical  analysis  methods  for  the  bending,  buckling,  and  vibration 
of  laminated  plating,  and  (c)  finite-element  techniques  for  either 
individual  laminated  components  or  more  complex  hybrid  structures. 
Finite-element  approaches  are  discussed  for  investigating  both  the 
in-plane  stresses  in  the  individual  layers  of  a  laminate  and  the 
interlaminar  shear  and  normal  stresses  which  frequently  govern 
near  discontinuities  in  laminate  cross  sections.  The  use  of  these 
three  analysis  techniques  necessitates  that  laminated  composite 
components  be  analyzed  as  nonlaminated  components  having  equivalent 
stiffness  during  in-plane  axial  and  shear  and  out-of-plane  bending 
responses.  Three  kinds  of  equivalent  elastic  properties  for 
laminates  are  defined,  and  methods  are  presented  for  calculating 
each  of  them.  Since  these  computational  procedures  may  be  quite 
lengthy  for  laminates  with  more  than  a  few  layers,  a  computer 
program  is  referenced  which  automates  the  procedures.  Lastly, 
the  effectiveness  of  different  methods  for  laminate  stiffness  and 
stress  analysis  was  investigated  by  comparing  analytical  and  test 
results  for  a  composite  box  beam  having  graphite  epoxy  flanges. 

Good  agreement  between  analytical  and  test  results  provides 
encouragement  for  the  further  application  of  these  methods. 


ADMINISTRATIVE  INFORMATION 

The  research  work  published  in  this  report  was  sponsored  during  fiscal  years 
1977-1980  by  the  Naval  Sea  Systems  Command  (Codes  05R  and  32R)  and  performed  at 
DTNSRDC  primarily  under  Work  Unit  1730-610. 


INTRODUCTION 

Composite  materials  have  been  employed  in  aerospace  structures  for  many  years 
and  are  now  receiving  increased  consideration  for  application  to  ship  structural 
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Figure  2  -  Composite  Foil  Flap  Component 


In  classifying  composite  materials,  three  commonly  accepted  types  are  recog¬ 
nized:  laminated  composites,  fibrous  composites,  and  particulate  composites.  Thi 

report  is  concerned  with  the  first  type,  the  laminated  composites,  which  are  norma ' 
fabricated  by  stacking  and  orienting  layers  of  fiber-imbedded  cloth  or  tape  in  pro 
ferred  directions  to  achieve  desired  strength  and  stiffness  properties.  An  examp  1 
of  a  four-layer  laminate  is  shown  in  Figure  3. 

As  illustrated  in  Figure  4,  existing  and  potential  applications  of  laminated 
composites  in  ship  structural  components  include:  (a)  laminated  beams  and  columns 
(b)  laminated  plating  used  as  the  flanges  of  beams  or  columns,  (c)  laminated  plate 
and  laminated  plating  for  (d)  panels  and  grillages,  and  (e)  the  skins  of  box-type 


structures . 


Numerous  methods,  carried  over  from  the  analysis  of  metallic  structures,  are 
available  for  analyzing  laminated  composites  in  the  above  applications.  These 
methods  include:  strength-of -mater ials  techniques  and  formulas;  classical  methods 
for  the  bending,  buckling,  and  vibration  of  laminated  plates;  and  the  finite-element 
method  for  the  analysis  of  more  complex  hybrid  structures  (see  Table  1) .  The  exten¬ 
sion  of  these  analysis  methods  from  structures  made  of  isotropic  and  homogeneous 
materials,  for  example,  metallic,  to  ones  involving  laminated  composites  necessitates 
that  the  laminated  composite  materials  be  replaced  by  "equivalent"  homogeneous  aniso¬ 
tropic  materials  having  the  same  stiffness  properties  during  extensional  and  bending 
deformations  (see  Figure  5).  The  required  analyses  are  then  performed  on  the  new 
structure  where  the  laminated  structural  components  are  assumed  to  be  made  of  this 
"equivalent"  material. 


LAMINATED  AND  EQUIVALENT  (NON LAMINATED)  BEAM  OR  COLUMN 
COLUMN 

(a) 


NOTE:  ABOVE  PROPORTIONS  FOR  Ef  <'  Ew 


BEAM  WITH  LAMINATED  FLANGES  AND  EQUIVALENT 
MODELS 

(bl 


Figure  5  -  Illustration  of  Modeling  Procedures  for  Laminated  Beams  and  Columns 
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The  first  part  of  this  report  will  present  procedures  ’  for  solving  for  the 
equivalent  stiffness  properties  based  on  the  principles  of  laminated  plate  theory.’*’ 
The  elastic  properties  needed  in  the  structural  analysis  of  laminates  are  of  three 
kinds  (Table  1):  (a)  engineering  constants  for  conventional  strength-of-mater ials- 

type  calculations,  (b)  extensional  and  bending  stiffnesses  for  classical  bending, 
buckling,  and  vibrational  analyses  of  laminated  plates,  and  (c)  reduced  stiffnesses 
for  finite-element  analyses  of  structural  components  involving  laminates.  The 
necessity  for  three  different  representations  of  laminate  stiffness,  and  the  distinc 
tion  between  them,  will  be  brought  out.  later  in  the  report. 

The  task  of  characterizing  elastic  properties  of  laminated  composites  may  be 
performed  analytically,  experimentally,  or  both.  The  use  of  laboratory  component 
tests  to  establish  the  elastic  properties  of  a  laminate  for  many  different  lamina 
layups  and  fiber  orientations  is  costly,  so  motivation  exists  for  having  an  experi¬ 
mentally  valid  analytical  tool  with  the  same  capability.  One  of  several  available 
computer  programs  having  this  analytical  capability  is  discussed  later  in  the  report 
This  and  the  other  programs  require  as  input  data  the  unidirectional  engineering 
constants  (parallel  and  perpendicular  to  the  fibers)  of  the  individual  layers  or 
laminas  which  make  up  a  laminate.  Although  these  unidirectional  engineering  con¬ 
stants  may  be  determined  analytically,  it  is  common  practice  to  experimentally 
measure  the  constants  using  coupon  tests. 


TABLE  1  -  EQUIVALENT  STIFFNESS  PROPERTIES  AND  METHODS 
OF  ANALYSIS  FOR  LAMINATED  COMPOSITES 


Methods  of  Analysis 

Types  of  Stiffness 

Proper t ies 

Strength  of  Materials 
(Formulas,  etc.) 

Engineering  Constants — 

E  ,  E  ,  v  ,  v  ,  G 
x  y  xy  yx  xv 

Classical  Plate  Solution 

Extensional  and  Bending 

Stiffnesses — A..,  B..,  D,. 

ij  ij  ij 

Finite  Element 

Reduced  Stiffnesses  Q*^ 

Q* .  ,  .-Membrane  Response 
ij(m)  p 

Q*  -Biding  Response 

Following  a  discussion  of  the  procedures  for  determining  the  elastic  properties 
of  laminates,  the  use  of  these  properties  in  implementing  the  analytical  methods  re¬ 
ferred  to  earlier  will  be  discussed.  In  order  to  assess  the  effectiveness  of  the 
strength-of-material  and  finite-element  methods  for  laminate  applications,  numerical 
and  test  results  are  reported  and  compared  for  a  composite  box  beam  simulating  the 
forward  foil  of  the  PCH-1  hydrofoil. 

EQUIVALENT  ELASTIC  PROPERTIES  OF  LAMINATES 

This  section  describes  the  computational  procedures  for  determining  the  equiv¬ 
alent  elastic  properties  of  laminates  which  are  needed  in  analyzing  structures  em¬ 
ploying  laminated  composites.  It  will  be  seen  that  these  various  kinds  of  elastic 
properties  are  interrelated  and  are  natural  to  the  particular  type  of  analysis  being 
performed.  Engineering  constants  are  common  to  design  and  analysis  methods  based  on 
strength-of-materials  techniques;  extensional  and  bending  stiffnesses  are  natural  to 
classical  methods  of  analyses,  including  finite-difference  techniques;  reduced  stiff¬ 
nesses  are  used  for  finite-element  analyses.  Because  the  engineering  constants  and 
reduced  stiffnesses  may  be  derived  from  the  extensional  and  bending  stiffnesses,  the 
latter  are  now  discussed. 

EXTENSIONAL  AND  BENDING  STIFFNESSES 

Mathematical  expressions  for  the  extensional  and  bending  stiffnesses  of  a 
laminate  may  be  derived  using  the  theory  of  laminated  plates,  or  lamination  theory 
as  it  is  called.^’*’  This  theory  is  empl  ved  to  not  only  generate  equivalent  stiff¬ 
ness  parameters  for  laminates,  but  also  to  form  the  basis  for  laminate  stress  anal¬ 
yses.  Lamination  theory  is  based  on  the  following  assumptions:  (a)  a  perfect  bond 
exists  between  the  layers  of  a  laminate  (no  slippage  between  layers) ;  (b)  the  lami¬ 
nated  plate  is  thin,  i.e.,  thickness  is  small  compared  to  the  lateral  dimensions  of 
the  plate;  and  (c)  normals  to  the  layers  of  a  laminate  remain  perpendicular  to  the 
layers  (shearing  strains  Yxz  =  YyZ  =  0)  and  do  not  change  length  (ez=0)  during 
laminate  bending.  The  coordinate  system  for  defining  ez,  Yxz»  and  YyZ  has  the 
z  axis  perpendicular  to  the  laminate  and  the  x  and  y  axes  lying  in  the  plane  of  the 
laminate  (see  Figure  6).  It  should  be  noted  that  assumptions  (b)  and  (c)  above 
correspond  to  the  Kirchoff  theory  for  plates  and  the  Kirchof f-Love  theories  for 
shells . 
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Figure  6  -  Relationship  between  (a)  In-Plane  Forces  and  Strains  and 
(b)  Moments  and  Curvatures  in  Flat  Symmetric  Laminates 


The  extensions  1  stiffnesses  A ^ .  of  a  laminate  relate  the  in— plane  forces 


■s'y >  and  ^’xy  per  unit  width  to  the  in-plano  strains  • 


>  and  o',  the 

v  XV 

middle  surface  (see  Figure  6a).  Similarly,  the  bending  stiffnesses  D..  relate  f, 

out-of-plane  bending  moments  M  and  M  and  twisting  moments  M  per  unit  width  to 

v  xv 


the  resulting  out-of-plane  curvatures  k  ,  k  ,  and  k  of  the  middle  surface  (see 

x  y  xy  ' 

Figure  6b) .  In  laminates  where  the  layup  of  layers  or  laminas  about  the  middle 
surface  is  not  symmetric,  either  because  of  fiber  orientation  and/or  material  pro 
erties,  a  coupling  between  extension  (including  shear)  and  bending  (including  twi 
ing)  occurs.  In  other  words,  the  application  of  in-plane  forces  not  only  produce 


y 


t 


in-plane  strains  but  also  causes  out-of-plane  bending  and  associated  curvatures. 
Likewise,  any  bending  moments  applied  to  the  laminate  will  produce  in-plane  strains. 
The  stiffness  parameters  which  relate  in-plane  forces  to  out-of-plane  curvatures  and 
which  relate  moments  to  in-plane  strains  are  referred  to  as  coupling  stiffnesses  B^. 
Mathematically,  these  relationships  are  expressed  in  matrix  form  by 


N 

X 

A11 

A12 

A16 

B11 

B12 

B16 

e 

X 

N 

y 

A12 

A22 

A26 

B12 

B22 

B26 

e 

y 

N 

xy 

A16 

A26 

A66 

B16 

B26 

B66 

Yxy 

M 

X 

B11 

B12 

B16 

D11 

°12 

D16 

k 

X 

M 

y 

B12 

B22 

B26 

°12 

D22 

D26 

k 

y 

M 

xy 

B16 

B26 

B66 

D16 

D26 

°66 

k 

xy 

(1) 


or,  in  condensed  form,  by 


(2) 


The  extensional,  coupling,  and  bending  stiffnesses  respectively,  and 

D„  in  Equation  (2)  are  expressed  and  determined  as  follows: 


-h/2 

N 

AU  =  J 

dz  =^’(Qij)k  (vVl* 

-h/2 

k=l 

f  h/2 

N 

BiJ  =  J 

ZdZ  =  2  (Zk_Zk-l) 

-h/2 

k=l 

h/2 

N 

=  J  *ij 

z2dz  ^(Qij)k  (zj-z IS 

-h/2 

k=l 
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where  the  Q  parameters  are  referred  to  as  reduced  stiffnesses  and  constitute  the 
elements  of  the  two-dimensional,  orthotropic  stress-strain  relations  for  the  individ¬ 
ual  layers  of  a  laminate.  These  stress-strain  relations  are  given  by 

(°x  \  **11  **12  **16  Cx  I 

°y  1  =  *12  *22  *26  «  Ey  >  (4) 

Txyj  _  **16  **26  **66_  ^xy / 

which  are  defined  with  respect  to  the  direction  of  the  loads  on  the  laminate.  The 

appearance  of  the  stiffness  elements  Q^,  ^26’  anc*  **66’  ■'-nstea^  ^13’  ^23’  an ^ 

follows  from  the  fact  that  Equations  (4)  are  obtained  by  reducing  the  6  >  6 

stress-strain  matrix  for  3-dimensional  behavior  (see  Reference  5).  Appendices  A 

and  B  and  Reference  5  discuss  the  detailed  procedure  for  determining  the  extensional 

and  bending  stiffnesses  A.,  and  D..  (B..  elements  are  zero  because  the  laminate  con- 

lj  ij  ij 

sidered  is  symmetric).  The  basic  steps  involved  are  summarized  as  follows: 

1.  Determine  the  reduced  stiffness  Q  for  each  of  the  laminas  (see  Appendix  A, 
Section  A.l)  using  available  unidirectional  elastic  moduli  E^,  E^,  ;21’  anc^  ''"l'1 

for  each  of  the  laminas. 

2.  Transform  these  Q...  defined  with  respect  to  the  so-called  material  direc- 

1J 

t ions  (parallel  and  normal  to  the  fibers),  to  the  directions  of  the  laminate  axes” 
x  and  y,  resulting  in  Q  (see  Appendix  A,  Section  A. 2). 

3.  Substitute  the  into  Equations  (3)  and  solve  for  the  A„  and  D^_.  (see 

Appendix  B)  . 

The  steps  iust  outlined  can  be  straightforwardly,  but  tediously,  used  to 

compute  the  required  stiffnesses  A^.  .  ,  B..,  and  D,.  for  any  given  laminate.  lhc 

i.l  1 J 

calculation  procedure  is  demonstrated  in  Appendix  C  for  a  graphite-epoxy  laminated 
skin  (symmetric  with  B  =  0)  of  the  composite  box  beam  in  Figure  1. 

At  this  point,  it  is  worthwhile  discussing  the  effect  of  such  laminate  char¬ 
acter  ist  ics  as  symmetry,  fiber  orientation,  etc.,  on  the  magnitudes  of  the  stillness 

coefficients  A..,  B..,  and  D..  in  Equation  (2).  Laminates  for  which  the  A,  B,  and  P 
ij  ij  1] 


*The  x  and  y  directions  for  a  laminate  are  usually  chosen  as  follows  tor  .1 
laminated  beam  and  column:  x  along  the  member  and  v  normal  to  the  member  in  the 
plane  of  the  laminate;  for  a  laminated  plate:  x  and  y  are  parallel  to  the  long  and 
short  sides  of  plate,  respectively;  for  laminates  which  are  components  of  2D  and.  o' 
hybrid  structures:  any  conveniently  oriented  x  and  v  axes  (z  too)  may  he  used. 


matrices  in  Equation  (2)  are  completely  full  (i.e.,  no  zero  elements)  represent  the 
most  general  case  encountered  (nonsymmetric  laminates  with  multiple  anisotropic 
layers) .  In  practice,  a  structural  designer  frequently  tries  to  avoid  this  general 
laminate  case  in  order  to  (a)  minimize  analytical  complexity,  (b)  eliminate  modes 
of  deformation  which  could  adversely  affect  the  load  capacity  of  the  laminated 
member  (e.g.,  out-of-plane  bending  and/or  twisting  deformations  associated  with  the 
coupling  stiffnesses  B^  will  reduce  the  load  capacity  of  a  laminated  plate  under 
in-plane  loads),  and  (c)  prevent  twisting  of  laminates  when  subject  to  temperature 
changes  and  in-plane  boundary  restraint.  Instead,  the  designer  often  tries  to 
select  or  buildup  a  laminate  having  symmetry  about  the  midplane  so  that  the  coupling 
coefficients  B_^  in  Equation  (2)  are  all  equal  to  zero.  Figure  7  shows  how  this  is 
achieved  for  cross-ply  and  angle-ply  laminates.  Note  that  symmetry  requires  that 
the  laminate  layers,  which  are  shown  reflected  (as  in  a  mirror)  about  the  middle 
plane,  must  have  the  same  thickness,  elastic  properties,  and  angles  of  orientation 
relative  to  the  x  and  y  axes. 

Beyond  the  simplification  in  Equations  (2)  due  to  symmetry  (B_  =  0),  simpli¬ 
fications  in  the  A  and  D  matrices  also  occur  if  additional  conditions  are  imposed 
on  the  arrangement  of  layers  within  the  laminate.  Figure  7  indicates  two  classes  of 
symmetric  laminates  which  may  be  somewhat  easily  discussed:  the  symmetric  cross-ply 
and  the  symmetric  angle-ply  laminate. 

For  the  case  of  symmetric  cross-ply  laminates,  whether  thay  are  regular  (e.g., 
90°,  0°,  90°)  or  irregular*  (e.g.,  90°,  0°,  0°,  90°)  in  layup,  the  coupling  coeffi¬ 
cients  A16,  A26,  D16,  and  D26  are  all  equal  to  zero.  Coupling  coefficients  A16  and 
A26  relate  extensional  forces  and  to  shear  strain  y ^  (and  shear  force  to 

extensional  strains  € ^  and  £  ) .  Similarly,  coupling  coefficients  and  D2&  relate 

bending  moments  M  and  M  to  twisting  curvature  k  ,  and  twisting  moment  M  to 
°  x  y  xy  xy 

bending  curvatures  k^  and  k^.  The  zeroing  out  of  A.^,  A^,  D^>  an<i  1*26  ^or  sym~ 
metric  cross-ply  laminates  follows  directly  from  the  fact  that  the  reduced  stiff¬ 
nesses  Q, ,  and  are  zero  for  all  of  the  layers  of  the  cross-ply  laminate. 

^16  26 

Turning  next  to  the  symmetric  angle-ply  laminates  in  Figure  7,  the  picture  is 

more  complicated.  For  the  regular  angle-ply  laminates  (e.g. ,  +a, -a,  +a) ,  the  A^, 

A„.,  D,,,  and  D„.  take  on  small  (nonzero)  values.  However,  for  the  irregular* 

26  16  26 


*Irregular  here  denotes  laminates  where  the  orientations  of  successive  layers 
do  not  simply  alternate  signs  or  directions  as  in  regular  laminates. 


12 


EXAMPLE-REGULAR  <90  ,  0  ,  90  ) 


COUPLING  COEFFICIENTS 


Bij  0  (DUE  TO  SYMMETRY) 

A16  A26  =  0 

°16  °26  ~  0 

FOR  REGULAR  OR  IRREGULAR 
CASES  (ANY  NUMBER  OF  LAYERS) 


EXAMPLE-REGULAR  (+a,  a,  +a)  OTHER  EXAMPLES 


REGULAR  LAYUPS 

A16'  A26'  °16'  °26 - ►TAKE  ON  SMALL 

VALUES-HR)  ABOVE 

IRREGULAR  LAYUPS 

A16'  A26'  D16'  °26 - ►TAKE  ON  SMALL 

VALUES— HI)  ABOVE 
A16'A26'D16'D26  =  0  «)'  ABOVE 


Figure  7  -  Coupling  Coefficients  for  (a)  Symmetric  Cross-Plv  (Specially 
Orthotropic)  and  (h)  Symmetric  Angle-Ply  Laminates 


symmetric  angle-ply  laminates,  the  coupling  coefficients  may  either  be  zero  or  non¬ 
zero  in  value  depending  upon  the  laminate  layup  arrangement  (orientations  of  the 
successive  layers)  as  pointed  out  in  Figure  7.  Lastly,  it  should  be  noted  that  a 
symmetric  laminate  put  together  from  a  combination  of  cross-plys  and  angle-plys  mav 
also  be  shown  to  lead  to  zero  on  nonzero  values  for  the  coupling  coefficients  depend¬ 
ing  on  the  arrangement  of  layers  within  the  laminate. 
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Recall  for  the  regular  symmetric  angle-ply  laminates,  that  the  coupling  coeffi¬ 
cients  A^,  A 2^ ,  D^,  and  D ^  took  on  small  nonzero  values.  In  these  instances  the 
analyst  should  be  cautioned  to  not  set  the  off-diagonal  terms  equal  to  zero  simply 
because  they  happen  to  be  small.  An  indication  of  the  effect  of  setting  the  off- 
diagonal  terms  and  equal  to  zero,  i.e.,  treating  the  angle-ply  laminated 
plate  as  specially  orthotropic,  is  provided  by  two  examples  from  Jones. ^  The  first 

example  is  a  uniformly  loaded  square  plate  for  which  D-./D,.  =  1,  (D, -  +  2D,,)/D11  = 

11  11  1 1  ob  11 

1.5,  and  D-^/D^  =  ^26^°11  =  For  cbis  problem,  setting  the  twist  coupling 

coefficients  and  equal  to  zero  results  in  calculated  plate  deflections  which 

are  24  percent  less  than  the  analytical  predictions  where  these  coefficients  have  not 

been  set  equal  to  zero.  The  second  example  in  Jones'*  has  to  do  with  the  buckling  of 

a  rectangular  laminated  plate  made  up  of  twenty  boron-epoxy  layers  having  alternating 

+ and  —angles  of  orientation  with  respect  to  the  overall  plate  axes.  It  is  pointed 

out  for  this  example  in  Jones  that  the  specially  orthotropic  approximation  leads  to 

an  overestimate  of  the  plate  buckling  load.  Additional  angle-ply  laminates  need  to 

be  looked  at  before  it  can  be  determined  just  how  small  the  D, and  D„.  have  to  be 

lb  lb 

before  they  can  be  safely  set  equal  to  zero  in  laminate  deflection  and  buckling 
analyses. 


ENGINEERING  CONSTANTS 

At  the  preliminary  design  stage,  many  aspects  of  structural  behavior  for  lami¬ 
nated  composite  members  (Figure  1)  may  be  initially  evaluated  using  strength-of- 
materials  techniques  once  the  necessary  engineering  constants  for  the  laminates  are 
available.  These  engineering  constants  may  be  calculated  from  the  extensional 
stiffnesses  A^.  after  their  values  have  been  determined  using  the  procedure  described 
previously  and  found  in  Appendices  A  and  B.  Methods  will  be  discussed  below  for 
calculating  the  following  engineering  constants  for  laminated  composites: 

E 

x 

v 

x 
G 

> 

n„ 

y>*y 

H  _  and  n  -  Coefficients  of  mutual  influence  of  the  second  kind 


1  and  E 

-  Young's  moduli 

x  y 

>  and  v 

-  Poisson's  ratios 

xy 

yx 

-  Shear  modulus 

’xy 

1  and 

r|  -  Coefficients 

x,xy 

y.xy 

1  and 

H  -  Coefficients 

xy  ,x 

xy,y 
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The  first  three  types  of  engineering  constants  above,  i.e.,  Young's  modulus, 
Poisson's  ratio,  and  shear  modulus,  are  already  familiar  to  those  individuals  in¬ 
volved  with  the  design  and  analysis  of  metallic  structures.  The  "coefficients  of 
mutual  influence,"  however,  are  most  probably  not  familiar.  Therefore,  stated 
simply,  the  coefficients  of  mutual  influence  of  the  first  kind  relate  a  normal  strain 

in  the  x  direction  to  a  shearing  strain  in  the  xy  plane  (p  )  or  relate  a  normal 

x ,  xy 

strain  in  the  y  direction  to  a  shearing  strain  in  the  xy  plane  (p  ).  Similarly, 

„  y,*y 

the  coefficients  of  mutual  influence  of  the  second  kind  relate  a  shearing  strain  in 

the  xy  plane  to  a  normal  strain  in  the  x  direction  (p  )  or  to  a  normal  strain  in 

xy,x 

the  y  direction  (p  ).  As  an  example,  p  is  expressed  bv 

xy,y  K  xy , x 


P„  =  Y  /c 

xy,x  xy  x 


where  Y  is  the  shearing  strain  in  the  xy  plane  associated  with  a  normal  strain 
in  the  x  direction  (usually  produced  by  a  normal  stress  0  )  . 

Methods  of  Calculation 

Several  computational  procedures  may  be  used  to  establish  the  engineering  con¬ 
stants  for  laminated  composite  materials.  One  procedure,  not  discussed  here,  is 

3  4 

contained  in  a  paper  by  Greszczuk.  A  second  procedure,  given  in  a  paper  by  Smith, 

is  probably  more  suitable  for  computer  implementation  than  the  first  and  is 

described  below. 

I'nder  in-plane  load  conditions,  the  in-plane  forces  on  the  laminated  compon¬ 
ents  in  structural  applications  (a)  through  (e)  of  Figure  4  are  related  to  the 
in-plane  strains  by 
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or,  in  matrix  form,  by 


{N}  =  [A]  (ej 


of  {N} 
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where  the  A.,  are  the  elements  of  the  matrix  [A]  which  is  the  inverse  of  [A]. 

A.,  are  found  from 
ij 


The 
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where 


A  An(A22  A66  ~  A26)  A12(A12  A66  A16  A26) 


+  A16(A12  A26  “  A16  A22) 


Once  the  A^  above  are  known,  these  quantities  may  be  used  to  compute  the  engineering 
constants.  Expressions  relating  the  engineering  constants  and  the  A^  are  obtained 
using  Equation  (6).  The  procedure  for  deriving  expressions  for  E^  and  ^  is  now 
shown.  (Similar  procedures  may  be  used  to  derive  expressions  for  the  other  engineer¬ 
ing  constants.) 
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Deriving  E^. 


Setting  =  N  =  0  in  the  first  of  Equations  (6),  we  have 

£  "A  N 
x  11  x 


e  -  o  / E  =  (N  /h)/E 

XXX  X  X 

Therefore,  equating  both  expressions  for  e  above  and  then  solving  for  E  , 

X  X 


Ex  ■  1/(Allh> 


where  h  is  the  full  thickness  of  the  laminate. 


Deriving  n 

_  *y,x 

Setting  =  N  =  0  in  the  third  of  Equations  (6), 

Y  =  A  ,  N 
xy  16  x 

Substituting  =  o^h  =  E^  e^h  into  the  above  equation, 

Y  =  A.  ,  E  eh 
xy  16  x  x 

Now,  as  defined  earlier, 

n  =  y  /e 
xy,x  xy  x 


A,  ,  E  eh 
16  x  x 


-  A.  ,  E  h 
16  x 


Therefore,  after  substituting  Equation  (8)  above  for  E  , 

n  =  a, ./a 

xy,x  16  11 


Having  shown  the  procedures  for  deriving  two  of  the  engineering  constants,  the 
final  expressions  for  all  of  the  engineering  constants  (including  those  just  derived) 
are  now  summarized: 
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■  1/<Allh> 


Ey  -  l/(*22h) 


V  -  -  e  /e  =  -  A.  ,,/A., 
xy  y  x  12  11 


v  =  -  c  /e  =  -  A10/A_0 
yx  x  y  12  22 


C*y  -  1/(A66h) 


(9) 


—  —  k 

n  =  a,, /a., 

xy,x  16  11 


n  =  a0,/a0„ 

xy,y  26  22 


n  =  a, ,/a,, 

x,xy  16  66 


H  =  A,,, /A,, 
y,xy  26  66 


A  simplification  of  the  process  for  calculating  engineering  constants  is 
possible  for  the  regular  symmetric  cross-ply  laminates  in  Figure  7. 

For  this  case,  since  A^  =  A^  =  0,  we  have,  using  Equations  (7)  and  (9), 


A  ~  A11  A22  A66  "  A12  A66 


and 


Therefore, 


A11  "  (A22  A66  A26^A  "  A22  A66/(A11  A22  A66  '  A12  A66) 


Ex  =  l/(A11h)  =  (An  -  A12/A22)/h 


(10) 


*References  4  and  5  only  included  expressions  for  the  coefficients  of  the 

second  kind  n  and  n  •  It  would  appear  that  these  coefficients  are  more 
xy , x  xy , y 

frequently  used  than  the  coefficients  of  the  first  kind  n  and  n 

7  x,xy  y,xy 
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Similarly,  it  may  be  shown  that 


and 


E 

y 


^A22  "  A12^All^h 


v 

xy 


A12^A22 


v 

yx 


A12/All 


G 

xy 


A66/h 


(11) 


n  =  n 

xy ,  x  xy , y 


n  =  n 

x,xy  y ,xy 


0 


It  should  again  be  emphasized  that  these  formulas  only  apply  for  regular  symmetric 
cross-ply  and  not  angle-ply  laminas.  For  angle  ply,  Equations  (8)  and  (9)  in  terms 
of  A. .  must  be  used. 


Experimental  Verification 

Using  the  two  methods  of  calculation  referred  to  in  the  previous  section,  engi¬ 
neering  constants  have  been  determined  for  a  laminated  composite  consisting  of  53 
plies  of  T300  graphite  fibers  oriented  0°  to  the  loading  direction  and  32  plies  of 
GY70  graphite  fibers  oriented  Z 45° .  Based  on  Method  I  (Smith's  paper  ),  engineering 

constants  are  calculated  in  Appendix  C  and  are  reported  in  the  first  column  of  Table 

3 

2.  Engineering  constants  based  on  Method  II  (Greszczuk's  work  )  were  taken  from  the 
theoretical  curves  of  Reference  1  (reproduced  in  Figures  8-12)  by  reading  off  the 
data  points  for  a  T300  and  GY70  hybrid  composition  of  55.3  percent  and  44.7  percent, 
respectively.  The  theoretically  predicted  engineering  constants  in  Table  2  are  seen 
to  be  in  close  agreement  with  the  experimentally  measured  values'*-  for  the  hybrid 
laminate.  It  is  concluded  that  either  Method  I  or  II  is  sufficiently  accurate  for 
calculating  the  engineering  constants  for  composit  laminates  under  in-plane  force 
and  strain  conditions.  Since  both  methods  are  mathematically  equivalent,  the  very 
slight  difference  in  numerical  values  is  most  probably  associated  with  (a)  interpo¬ 
lating  the  theoretical  values  from  the  theoretical  curves  in  Figures  8-12  for 
specified  percentages  of  T300  and  GY70  and  (b)  plotting  up  the  theoretical  curves  in 
Figures  8-12  in  the  first  place. 
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TABLE  2  -  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  STIFFNESS 
PROPERTIES  FOR  BOX  BEAM  HYBRID  COMPOSITE 
(THORNEL  T300  AND  CELION  GY70) 


Engineering 

Theoretical 

Theoretical 

3 

Constant 

(Method  I)*4 

(Method  II)3 

Experimental 

E 

X 

12.67 

12.6 

13.04 

E 

4 . 64 

4.6 

4.17 

y 

V 

0.804 

0.800 

0.766 

xy 

V 

yx 

0.294 

0.295 

0.244 

G 

5.26 

5.2 

5.09 

xy 

COMPUTER  PROGRAM 

Since  the  procedure  for  the  calculation  of  extensional,  bending,  and  coupling 
stiffnesses,  as  well  as  the  engineering  constant,  is  a  rather  tedious  one,  especially 
for  a  composite  laminate  involving  many  layers,  computer  programs  have  been  written 
to  automate  the  process.  One  such  program^  is  called  SQ5  and  is  available  through 
the  Aerospace  Structures  Information  and  Analysis  Center,  Wright-Patterson  Air  Force 
Base,  Ohio.  This  program  calculates  the  laminate  stiffness  properties  mentioned 
above  and  performs  a  stress  analysis  for  a  given  set  of  applied  in-plane  forces  and 
moments. 

5  6 

The  theoretical  basis  for  SQ5  is  plate  lamination  theory.  ’  First,  for  a 
given  laminate,  the  elements  of  matrix  Equation  (1)  are  generated.  These  equations 
are  then  solved  for  the  in-plane  strains  and  curvatures  in  terms  of  applied  forces 
and  moments  by  inverting  the  A  and  D  matrices  (B„  =  0).  The  strains  and  curvatures 
are  next  used  to  determine  the  strains  and  stresses  in  each  lamina.  Temperature- 
induced  stress  and  m  ment  resultants  may  also  be  calculated  and  then  added  to  the 
other  known  loads.  The  thermal  analysis  assumes  a  constant  temperature  throughout 
the  laminate  thickness.  The  program  also  includes  a  simplified  transverse  shear 
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♦  45  MATERIAL  (%) 


Figure  8  -  Young's  Modules  of  Laminate  as  Function  of  Percent  of  GY70 
(+45  Degrees)  and  T300  (0  Degree,  90  Degrees)  Material  (from  Reference  1) 


T 


±  45’  MATERIAL  (%) 

Figure  10  -  Poisson's  Ratio  of  Laminate  as  Function  of  Percent  of  GY70 
(+45  Degrees)  and  T300  (0  Degree,  90  Degrees)  Material  (from  Reference  1) 

23 
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±45°MATERIAL  (%) 

Figure  11  -  Poisson's  Ratio  V  of  Laminate  as  Function  of  Percent  of  GY70 

yx 

(+45  Degrees)  and  T300  (0  Degree,  90  Degrees)  Material  (from  Reference  1) 
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analysis  in  which  the  shear  stress  distribution  across  the  laminate  thickness  due  to 

known  shear  resultants  Q  and  Q  is  predicted.  It  also  has  the  capability  to  obtain 

x  y 

a  laminate  interaction  diagram  based  upon  the  maximum  strain  theory  of  failure. 

A  description  of  the  input  data  needed  to  use  the  program  SQ5,  excerpted  from 
Reference  7,  is  provided  in  Appendix  D. 

ANALYSIS  OF  LAMINATES 

STRENGTH-OF-MATERIALS  METHODS 

Strength-of-materials  methods  and  formulas  are  useful  in  varying  degrees  for  all 
of  the  structural  applications  of  Table  1  once  the  appropriate  equivalent  elastic 
properties  for  the  laminate  have  been  found.  However,  these  methods  are  probably 
most  effective  for  structural  applications  (a)  and  (b)  in  Figure  4  which  are  con¬ 
cerned  with  laminated  beams  and  columns.  In  application  (a),  the  cross  section  of 
the  member  is  considered  to  be  laminated  from  top  to  bottom  and  has  a  rectangular 
cross  section  (see  Figure  5).  In  application  (b) ,  only  the  flange  or  flanges  of  the 
member  are  of  laminated  construction  (Figure  5).  The  web  consists  of  a  homogeneous 
and  isotropic  metallic  or  nonmetallic  material. 

Laminated  Beams  and  Columns 

In  fabricating  beams  or  columns,  it  is  advantageous  to  stack  the  layers  so  as 
to  provide  symmetry  about  the  midplane  and,  thereby,  eliminate  coupling  between  ex- 
tensional  and  bending  responses,  as  discussed  earlier.  All  layers  should  probably 
be  oriented  with  the  fibers  parallel  to  the  length  of  the  beam.  In  some  engineering 
situations,  a  laminated  member  may  not  be  fabricated  from  scratch  but  may  be  ob¬ 
tained  simply  by  cutting  a  narrow  strip  from  a  laminated  plate  which  is  readily 
available  and  may  have  been  laid  up  for  other  purposes.  In  these  cases,  the  orien¬ 
tation  of  fibers  in  the  individual  layers  of  the  laminate  may  be  oriented  in  either 
a  cross-ply  or  angle-ply  arrangement,  or  a  combination  of  both. 
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The  governing  differential  equation  for  the  deflection  of  a  laminated  beam  or 
plate  strip*  (Figure  4a)  is  given  by 


E  1  ^  =  M 
b  2  x 
dx 


(12) 


where  y  is  beam  deflection,  is  bending  moment,  E^  is  the  engineering  constant  for 
an  equivalent  beam  having  the  same  bending  rigidity  E^I  as  the  actual  laminated  beam 
(see  Figure  5a),  and  I  is  cross-sectional  moment  of  inertia  (the  equivalent  beam  is 
assumed  to  have  the  same  cross-sectional  dimensions,  and  therefore  moment  of  inertia 
I,  as  the  laminated  beam). 

In  order  to  obtain  beam  deflections  by  integrating  Equation  (12)  and  satisfying 
boundary  conditions,  the  engineering  constant  E^  for  the  equivalent  beam  must  be 
determined.  To  do  this,  one  starts  with  the  moment-curvature  relations,  in  inverted 
form,  for  a  laminated  plate 


k 

x 


k 

y 


k 

xy 


Dln  M  +  D. _  M  +  D  ,  M 
11  x  12  y  16  xy 


D10  M  +  D._  M  +  D0,  M 
12  x  22  y  26  xy 


D.,  M  +  D„,  M  +  D,,  M 
16  x  26  y  66  xy 


(13) 


which  are  analogous  to  Equations  (6)  involving  in-plane  forces  and  strains.  Setting 
My  =  M  =  0  in  the  first  of  Equations  (13),  to  represent  one-dimensional  (beam) 
bend ing. 


or 


k 

x 


'11 


M 

x 


1 


k 

x 


M 

x 


(14) 


*The  plate  strip  here  is  assumed  to  be  in  a  state  of  plane  stress.  If  it  is 
desired  to  model  the  plate  strip  as  an  element  of  an  entire  plate,  it  is  necessarv  to 
assume  a  plane-strain  condition  for  the  sides  of  the  strip  and  to  start  with  the 
appropriate  differential  equation. 
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which  applies  for  a  beam  or  strip  of  unit  width.  For  a  beam  of  width  b,  we  have 


Du dx  x 


(15) 


Comparing  Equations  (15)  and  (12),  we  see  that  the  equivalent  engineering  constant 


is  obtained  from 


V  ■  b/DU 


(' '  & 


bh~ 


or 


(fc 


Eb  (  bh- 


3  i  .  b 


'll 


from  which 


Efc  -  12/ (Du  hJ) 


(16) 


It  should  be  noted  here  that  E^  associated  with  the  bending  response  of  a 

laminate  is  different  from  E^  in  Equation  (8)  associated  with  the  extensional 

response.  This  relates  to  the  fact  that  E^  is  dependent  upon  the  stacking  sequence, 

but  E  is  not. 
x 

The  bending  stiffness  D  in  Equation  (16)  is  needed  to  solve  for  E,  .  The 
_  -Li.  b 

stiffness  is  given  by 


D11  ■  <022  “66  -  D26)/4 


where 


A  ‘  D1I(D22  D66  -  D26>  '  D12<D12  °66  ‘  »16  “26' 


+  “l6<D12  D26  -  V  D22> 
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For  regular  symmetric  cross-ply  laminates  (Figure  7),  =  0. 

Therefore, 


11 


D11  -  D12/D22 


and 


Eb  = 


12  (D11  -  w 


(17) 


In  addition  to  determining  deflections  of  a  laminated  beam  using  Equations  (12) 
and  (16),  interest  also  exists  in  determining  stresses.  A  coarse  estimate  of  stress 
in  the  laminated  beam  is  obtained  from  the  formula 


o 

x 


M  c 
x 

I 


(18) 


where  M  is  the  bending  moment,  c  is  the  distance  from  the  neutral  surface  to  the 
fiber  location,  and  I  is  the  beam  cross-sectional  inertia. 

A  more  accurate  estimate  of  stress  in  any  particular  layer  of  the  beam  may  be 
obtained  as  follows.  Using  Equation  (18),  the  longitudinal  strain  in  any  layer  of 
the  beam  is  obtained  from 


£ 

X 


a  /E, 
x  b 


M  c 
x 


V 


The  lateral  strain  in  the  beam  is  found  using 


r.  =  -  v  e 

y  xy  x 


(19) 


(20) 


where  v  is  given  by  one  of  Equations  (9)  for  a  laminate  under  extension.*  Then, 
xy 

the  stress  in  any  layer  is  given  by 


*We  are  neglecting  the  fact  that 
extension  in  obtaining  v  .  In  order 

could  be  replaced  by  E^  from  Equation 


the  laminate  is  really  under  bending  and  not 
to  be  more  consistent,  E^  in  Equation  (19) 

(8). 
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(21) 


Q..  e  +  Q, »  e 
^11  x  12  y 


a  =  Q10  e  +  Q00  e 
y  12  x  22  y 


where  the  reduced  stiffnesses  Q-q«  an<^  ®22  are  ^oun<^  f°r  any  layer  or  lamina  of 

the  beam  using  the  theory  and  procedures  given  in  the  first  two  sections  of  Appendix 
A.  The  stress  0 x  will  probably  dominate  over  o  . 

Equation  (16)  for  is  also  needed  for  determining  the  buckling  load  of  a 

laminated  column  or  plate  strip.  The  critical  buckling  load  is  given  by  the 
Euler  formula 

P  =  tt2E,  I/L2  (22) 

cr  b 

for  a  pinned-end  laminated  column  of  length  L  and  cross-sectional  moment  of  inertia 


Members  with  Laminated  Flanges 

It  is  anticipated  that  laminated  members  may  involve,  in  some  instances,  a 
homogeneous  metallic  web  with  a  laminated  flange  as  illustrated  in  Figure  5b.  In 
addition  to  the  I  beam  shown,  the  hybrid  structural  member  may  also  take  the  form  of 
a  T  beam.  An  approximate  determination  of  the  flexural  rigidity  E^I  for  these  hybrid 
members  may  be  accomplished  using  a  two-step  process  as  shown  in  Figure  5b.  In  the 
first  step,  the  hybrid  I-member  with  laminated  flanges  is  converted  into  a  hybrid 
I-member  with  homogeneous  flanges  by  mathematically  replacing  the  laminated  flange 
material  with  a  nonlaminated  material  having  an  equivalent  engineering  constant  E^ . 
The  equivalent  constant  is  found  using  Equations  (7)  and  (8)  given  earlier.  The 
second  step  is  to  convert  the  hybrid  member  with  a  homogeneous  flange  material  to  a 
nonhybrid  member  where  the  flange  material  is  replaced  by  a  new  flange  of  width  b0, 
having  a  modulus  E^  identical  to  that  of  the  web  of  the  I  beam.  The  new  required 
flange  with  b_  is  found  from 


*See  a  strength-of-material  text  or  handbook  for  dealing  with  other  end  bound¬ 
ary  conditions. 


•  •  -*><r 


Once  the  initial  hybrid  beam  has  been  converted  into  the  final  equivalent  beam  just 
described,  the  moment  of  inertia  I ^  is  computed  for  this  final  cross  section.  The 
flexural  rigidity  E^^  *-s  then  substituted  for  E^I  in  Equation  (12)  for  calculating 
beam  deflections  and  in  Equations  (19)  as  the  first  step  in  determining  stresses  in 
the  individual  layers  of  the  laminated  flange.  Finally,  it  should  be  pointed  out 
that  the  analysis  procedure  for  a  structural  member  with  a  laminated  flange  intro¬ 
duces  an  approximation  which  is  not  present  in  the  laminated  beam  analysis  of 
Figure  5a.  The  approximation  is  associated  with  computing  the  equivalent  flange 
modulus  E^  above  by  assuming  that  the  strain  distribution  over  the  depth  of  the 
flange  is  uniform.  Actually,  the  strain  distribution  varies  slightly  due  to  bending 
of  the  overall  beam  cross  section. 

CLASSICAL  METHODS 

Laminated  plating  has  potential  applications  in  panels,  grillages,  and  three- 
dimensional  hybrid  structures,  as  illustrated  in  Figure  4.  For  analysis  purposes, 
the  plating  in  these  applications  may  often  be  modeled  as  a  single  plate,  as  shown 

in  Figure  13,  if  appropriate  edge  boundary  conditions  are  selected. 

8  5 

Timoshenko  and  Jones  discuss  classical  methods  for  analyzing  the  bending, 
buckling,  and  vibration  of  rectangular  orthotropic  plates  which  may  be  directly 
applied  to  regular  symmetric  cross-ply  laminates  (see  Figure  7) .  One  of  the  approxi¬ 
mate  solution  techniques  described  in  these  references  uses  Fourier  series  to  repre¬ 
sent  the  plate  loading  and  deflection  shape.  Formulas  are  presented  in  the  next 
section  which  were  derived  based  on  this  approach.  As  mentioned  above,  the  formulas 
apply  to  specially  orthotropic  laminates  involving  cross-ply  fibers  for  which  the 
stiffnesses  and  in  Equations  (l)-(3)  are  equal  to  zero.  The  solution  pro¬ 

cedure  for  regular  symmetric  angle-ply  laminates  is  more  complex  since  the  and 
D„,  stiffnesses  are  no  longer  zero  but  may  be  small  when  a  large  number  of  layers 
are  involved.  The  analysis  of  nonsymmetr ic  laminas  involving  nonzero  B_  coupling 
stiffnesses  between  extension  and  bending  introduces  even  further  complexity  into 
the  analysis. 
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Figure  13  -  Simply  Supported  Laminated  Plates  Subjected  to  (a)  Lateral 
Bending  and  (b)  Buckling  Under  In-Plane  Loads 

Lateral  Bending  of  Laminated  Plates 

Q 

Using  the  solution  procedure  of  Navier  which  is  described  by  Timoshenko  and 
Jones,  the  deflection  of  a  uniformly  loaded  plate  with  simply  supported  edges  is 
given  by 


lbp 


w(x,y)  =  — ^ 


Z  L 


1  .  /  mTTx  \  .  /  nTTy  \ 

—  sin  - 1  sin 

mn  '  a  '  '  n  / 


m-1,3,5,...  n-1,3,5,...  D„  ( +  2  (D,„  +  2D„)  (?)  +  f 


I 


where  is  the  uniform  load  per  unit  area, 

a  and  b  are  the  side  dimensions  of  the  plate, 

x  and  y  are  variable  distances  parallel  to  a  and  b  sides  of  plate, 
m  and  n  are  odd  integers,  and 

Dll’  D12’  D22  an<*  D66  are  in  Equation  (3)  and  computed  in  Appendix  B 

for  a  given  example  of  a  four-layer  laminate. 

The  strains  in  any  layer  of  the  laminate  at  a  distance  z  from  the  middle  surface 

may  be  found  by  substituting  Equation  (23)  for  w  into 


e 

x 


£  = 
y 


(24) 


and 


■xy 


-  z 


.2 

3  w 

3x9y 


Stresses  are  then  obtained  by  substituting  Equations  (24)  for  strain  into  Equations 
(4).  By  following  the  procedure  just  outlined,  formulas  have  been  generated  for 
computing  these  stresses.  In  arriving  at  these  formulas,  the  series  in  Equation  (23) 
has  been  written  out  for  m  taking  on  values  of  1  and  3  and  for  n  taking  on  values  of 
1  and  3.  The  stresses  in  any  layer  are  then  given  by 


a  =  Q, ,  e  +  Q.  „  e 
x  X11  x  12  y 


and 


o  =  e  +  Q„.  e 
y  x12  x  22  y 


(23) 


where 


16Pnz  „  „ 

ex  =  — —  ((C1+C2)(7T/a)Z  +  (C3+CA)(37T/a)^} 

7T 


and 


16p0z 


e  =  •— r-  {(C^+C  )(TT/b)2  +  (C2+C  )  (3TT/b)2 } 

y  TT 
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and 


C1  =  L/{Dn(l/a)4  +  2(D12+D66)(l/a)2(l/b)2  +  D22(l/b)4} 

C2  =  (-l/3)/{Dn(l/a)4  +  2(D12+D66)(l/a)2(3/b)2  +  D22(3/b)4} 
C3  =  (-l/3)/{Dn(3/a)4  +  2(D12+D66)(3/a)2(l/b)2  +  D22(l/b)4} 
C4  =  (  l/9)/{Du(3/a)4  +  2(D12+D66)(3/a)2(3/b)2  +  D^O/b)4} 


The  stiffnesses  Q^,  Q-^,  an<^  ^22  are  ca^cu^ate(^  f°r  eac^  layer  using  equations  in 
Appendix  A. 


Buckling  and  Vibration 

Equations  are  presented  in  this  section  for  calculating  the  buckling  loads  and 
vibration  frequencies  of  rectangular  laminated  plates  having  simply  supported 
boundaries  and  a  regular  symmetric  cross-ply  layup  of  fibers  (defined  as  specially 
orthotropic  plate,  see  Figure  7) . 

The  buckling  stress  0 ^  of  a  specially  orthotropic  plate  under  an  axial  loading 
5  8  9 

(Figure  13)  is  given  by  ’  ’ 


where  a,  b,  and  h  are  the  length,  width,  and  thickness  of  the  plate,  m  and  n  are 
integers  indicating  the  number  of  buckle  half  wavelengths  in  the  x  and  y  directions, 
and  ,  D^2,  and  D22  and  ate  defined  by  Equations  (l)-(3).  Rearranging  Equa¬ 
tion  (26)  and  recognizing  that  the  lowest  value  of  buckling  stress  will  be  associ¬ 
ated  with  n  =  1, 


a 

cr 


+  2(D12+2D66)  +  ° 


22 


(27) 
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The  number  m  of  half  waves  in  the  x  direction  is  determined  by  taking  the  deriv 
ative  do^/dm,  setting  it  equal  to  zero,  and  solving  for  m.  This  gives 


m 


■f 


(28) 


If  the  m  computed  using  Equation  (28)  happens  to  be  an  integer,  which  is  unlikely, 
this  value  of  m  gives  the  actual  number  of  half  waves  in  the  x  direction,  and  the 
corresponding  minimum  buckling  stress  results  from  substituting  Equation  (28)  into 
Equation  (29),  which  gives 


<0cr>»i„  -  + 


(I,i2+2D66>i  fr 

b  n 


(29) 


However,  if  the  m  computed  from  Equation  (28)  is  not  an  integer — the  most  frequent 
case — it  is  necessary  to  substitute  the  next  higher  and  lower  integer  values  into 
Equation  (27)  to  see  which  corresponds  to  the  lower  buckling  stress.  For  example, 
if  m  is  computed  to  be  1.6,  one  should  try  m  =  1  and  m  =  2  in  Equation  (26)  to  see 
which  m  values  results  in  the  lower  buckling  stress.  If  m  should  turn  out  to  be 
less  than  1,  say  m  =  0.6,  then  m  =  1  corresponds  to  the  lowest  buckling  mode  and 
the  corresponding  critical  stress  is  found  by  substituting  m  =  1  into  Equation  (27). 

Buckling  loads  for  specially  orthotropic  plates  under  other  loading  conditions 
(biaxial  and  shear)  as  well  as  other  boundary  conditions  (such  as  two  edges  ly 

supported  and  two  edges  clamped  or  elastically  restrained,  or  other  combinations) 
are  provided  in  Reference  9. 

In  addition  to  having  analytical  means  for  investigating  the  stability  of 
laminated  plates,  tools  are  also  needed  to  estimate  the  fundamental  frequencies  of 
vibration  of  laminated  plates.  Fundamental  plate  frequencies  are  necessary  to 
insure  that  these  frequencies  are  sufficiently  removed  from  exciting  frequencies 
due  to  machinery  or  other  cyclic  driving  sources  so  as  to  prevent  resonance  or 
near-resonance  conditions  which  may  result  in  excessive  deflections,  fatigue,  or 
other  structural  damage.  The  natural  frequencies  of  vibration  w  of  rectangular 
laminated  plates  with  simply  supported  boundaries  may  readily  be  computed  using 
the  equation 


35 


»n  (l)  *  (!)  (!)  +  ”22  (S 


where  p  Is  the  mass  density  of  the  plate,  a,  b,  and  D^'s  are  identical  to  those 
defined  for  laterally  loaded  plates  and  buckling,  and  ni  and  n  refer  to  the  different 
mode  shapes  of  vibration.  The  fundamental  frequency  of  vibration  is  obtained  by 
setting  m  and  n  equal  to  1.  Reference  6  gives  formulas  for  computing  the  natural 
frequencies  of  vibration  for  specially  orthotropic  laminated  plates  with  all  clamped 
boundaries  or  for  various  combinations  of  simply  supported  and  clamped. 


FINITE-ELEMENT  METHODS 

Analytical  methods  were  given  in  the  previous  sections  for  analyzing  several 
basic  laminated  components  in  hybrid  structures.  These  methods  isolate  a  laminated 
beam,  column,  or  plate  component  of  the  structure  and  introduce  the  appropriate 
boundary  and  loading  conditions  to  represent  the  interaction  of  the  isolated  member 
with  the  remainder  of  the  structure.  Since  the  results  oi  these  methods  (strength- 
of -materia Is  and  classical  analyses)  are  usually  of  an  approximate  nature,  the 
methods  are  most  suitable  for  application  in  the  early  stages  of  design.  In  order 
to  obtain  more  accurate  results  in  the  later  and  final  stages  of  design,  the  analyst 
frequently  turns  to  the  finite-element  method  and  the  computer  programs  which  imple¬ 
ment  this  method.  One  of  the  advantages  of  finite-element  analyses  is  that  these 
analyses  minimize  the  number  of  assumptions  as  to  the  loading  and  constraints  on  the 
various  components  of  the  structure. 

In  analyzing  hybrid  structures,  at  least  four  options  are  available  for  model¬ 
ing  an.i  treating  the  laminated  components  of  the  structure.  These  four  approaches 
are  illustrated  in  Figure  14  for  the  tapered  composLit  box  beam  of  Figure  1  which 
simulates  the  forward  roil  on  the  PCH-1  hydrofoil.  1  he  composite  box  beam  has 
graphite  epoxy  laminated  skins  with  internal  steel  spars  and  is  supported  and  loaded 
as  indicated  in  Figure  1. 

Approach  (l)  in  ‘•igure  14  entails  modeling  the  full  laminated  skin  thickness  in 
t he  box  beam  by  equivalent  two-dimensional  (2D)  piate  elements  having  anisotropic 
material  properties  which  are  equivalent  to  those  of  the  multilayered  laminate. 
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(2)  THREE-DIMENSIONAL  FINITE 
ELEMENTS  PROVIDE 
•IN-PLANE  STRESSES  (cT,  oy>  rxy) 

•INTERLAMINAR  AND  NORMAL 
STRESSES  (txz,  ryz,  ogi 


(3)  TWO-DIMENSIONAL  FINITE 

ELEMENTS  (yz  PLANE)  PROVIDE 
•IN-PLANE  STRESSES  (oy,  0x  =  0) 
•INTERLAMINAR  AND  NORMAL 
STRESSES  (rzyANDaz) 


(4)  THREE-DIMENSIONAL  FINITE 
ELEMENTS  PROVIDE 
•IN-PLANE  STRESSES  (oy  AND  oj 

•INTERLAMINAR  AND  NORMAL 
STRESSES  (rzy,  rzx,  AND  az> 


Figure  14  -  Finite-Element  Approaches  for  Laminated  Skins  in 

Composite  Box  Beam 
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This  modeling  procedure  was  employed  at  DTNSRDC*  in  analyzing  the  composite  box  beam. 

The  procedure  allows  a  determination  of  the  in-plane  normal  stresses  o  and  o  and 

x  y 

shear  stress  t  in  the  laminate,  in  an  average  or  equivalent  sense.  These  equiva¬ 
lent  laminate  stresses  may  be  used  to  determine  the  in-plane  stresses  in  the  individ¬ 
ual  layers  of  the  laminate  in  an  approximate  way  (see  Appendix  F) .  However,  it  is 
not  possible  to  determine  the  interlaminar  shear  stresses  and  t  by  adopting 
equivalent  2D  plate  elements  for  the  full  laminate.  The  latter  stresses  are  fre¬ 
quently  the  governing  ones,  particularly  near  free  edges  in  laminates  produced  by 
cutting  or  drilling  holes  during  fabrication. 

Approaches  (2)  through  (4)  in  Figure  14  facilitate  a  determination  of  interlami¬ 
nar  shear  and  normal  stresses.  In  applying  approach  (2)  in  Figure  14,  the  laminate 
thickness  is  discretized  by  stacking  a  suitable  number  of  3-dimensional  (3D)  solid 
elements  to  represent  the  full  thickness  of  the  laminate.  It  is  clearly  not  prac¬ 
tical  to  have  as  many  3D  elements  vertically  through  the  thickness  as  there  are 
actual  layers  since  there  are  frequently  a  large  number  of  layers  involved  (approxi¬ 
mately  80  layers  for  the  skin  of  the  box  beam).  Instead,  the  layers,  may  be  grouped 
into  as  many  layers  as  necessary  (four  layers  in  Figure  14)  to  determine  the  inter¬ 
laminar  shear  and  normal  stresses  at  desired  locations  between  layer  or  plies  of  the 
laminate.  This  second  modeling  approach  was  adopted  by  the  Virginia  Polytechnic 
Institute10  (VPI)  in  analyzing  the  composite  box  beam.  However,  using  this  approach, 
the  overall  box  beam  deflection  results  were  in  agreement  with  neither  the  results 
of  DTNSRDC's  f in it e- element  analyses  using  approach  (1)  nor  with  the  box  beam  test 
results.  It  is  suspected  that  the  major  reason  that  VPI's  analysis  did  not  agree 
with  the  others  was  their  use  of  finite  elements  having  very  large  aspect  ratios. 

In  addition,  the  distribution  of  interlaminar  shear  and  normal  stresses  along  the 
length  of  the  box  beam  also  exhibited  some  questionable  oscillations.  Therefore, 
further  work  is  apparently  still  needed  on  the  implementation  of  this  approach. 

Approaches  (3)  and  (4)  in  Figure  14  utilize  a  beam-type  model  of  the  box  beam. 
Because  of  this  modeling  approximation,  the  resulting  interlaminar  shear  and  normal 
stresses  are  not  as  accurate  as  they  would  be  with  approach  (2)  assuming  all  analyses 
are  carried  out  properly.  By  stacking  2D-type  elements  in  the  yz  plane  in  approach 
(3),  a  plane  stress  analysis  is  accomplished  for  which  0^  =  0,  as  indicated  in 


*Stein,  M.C.,  "A  Nastran  Analysis  of  a  Composite  Laminate  Box  Beam  for  Applica¬ 
tion  to  Navy  Hydrofoils,"  reported  informally  as  enclosure  (1)  to  DTNSRDC  ltr 
77-173-186  of  9  Dec  1977. 
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Figure  14.  Since  the  box  beam  skins  are,  to  some  degree,  restrained  in  the  trans¬ 
verse  direction  by  the  load  block  at  the  one  end  and  the  clamping  support  at  the 
other,  a  plane  strain  analysis  where  the  transverse  strains  in  the  skin  are  set 
equal  to  zero  (but  is  not  zero)  may  be  a  better  approximation  to  the  box  be¬ 
havior  when  using  a  beam-type  model.  This  plane  strain  analysis  may  be  accomplished 
by  stacking  3D  finite  elements  as  indicated  in  approach  (4)  of  Figure  14.  It  should 
be  noted  that  in  implementing  approaches  (3)  and  (4) ,  it  is  important  that  the 
portion  of  the  total  box  beam  load  which  is  carried  by  the  beam  strip  be  approxi¬ 
mately  determined  and  applied  to  tne  beam  model.  It  also  should  be  noted  that  both 
approaches  (2)  and  (4),  involving  3D  elements,  necessitate  that  five  additional 
engineering  constants  (over  those  required  with  the  2D  finite  elements  in  approaches 
(1)  and  (3))  must  be  experimentally  measured:  E  ,  V  ,  V  ,  G  ,  and  G  .  Finally, 
in  terms  of  the  cost  to  carry  out  the  above  approaches,  it  can  be  said  that  approach 
(3)  will  be  the  least  costly  and  (2)  the  most  costly.  Approaches  (1)  and  (4)  will  be 
somewhere  in  between  in  cost;  it  is  difficult  to  rank  them  with  respect  to  each  other 
since  the  cost  would  be  dependent  on  the  number  of  elements  used  in  each  case. 

Composite  Box  Beam  Analysis 

Using  the  computer  program  NASTRAN,11  a  finite-element  analysis  was  performed 
by  DTNSRDC  on  the  composite  box  beam  in  Figure  1  to  verify  the  results  of  a  strength- 
of-materials  analysis  and  to  provide  a  more  detailed  insight  into  the  beams  struc¬ 
tural  behavior.  The  NASTRAN  idealization  for  the  box  beam  is  shown  in  Figure  15. 

This  report  focuses  on  those  aspects  of  the  f inite— element  idealization  and  analysis 
concerned  with  the  laminated  skins  since  the  steel  spars  may  be  modeled  using  pro¬ 
cedures  already  familiar  to  the  average  NASTRAN  user. 

Modeling  of  Laminated  Skin.  The  skin  of  the  box  beam  is  a  reinforced  epoxy  laminate 
consisting  of  44.6%  GY70  at  +45°  and  55.4%  of  T300  at  0°,  with  a  constant  thickness 
of  0.5  in.,  except  for  regions  near  the  load  block  and  fixed  end.  The  modeling  of 
the  laminated  skin  was  accomplished  using  the  NASTRAN  plate  element  CQUAD1.  Accord¬ 
ing  to  the  NASTRAN  user's  manual,11  the  CQUADl  element  is  intended  for  application 
to  sandwich  plates  having  different  elastic  properties  in  bending,  membrane,  and 
shear.  No  mention  is  made  of  using  CQUADl  elements  for  laminated  plates  in  the 
manual.  However,  the  idea  occurred  at  the  beginning  of  the  composite  box  beam 
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CQUAD1  PLATE  MEMBRANE  AND 
BENDING  ELEMENTS 
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Figure  15  -  NASTRAN  Finite-Element  Model  for  Composite  Box  Beam 


analysis  that  the  CQUAD1  element  could  be  used  to  model  the  laminated  skin  as  an 
equivalent  (nonlaminated)  skin  having  the  same  thickness  and  same  anisotropic  elastic 
properties  under  membrane  and  bending  action.  (Transverse  shear  stiffness  was 
assumed  to  be  infinite.)  This  approach  was  therefore  followed  in  the  box  beam 
analysis . 


L 
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From  a  structural  point  of  view,  the  NASTRAN  membrane  plate  element  CQDMEM 
could  have  been  used  for  modeling  the  box  beam  skin  since  the  local  bending  of  the 
skin  is  minimal  due  to  the  configuration  and  loading.  Instead,  the  CQUAD1  element 
with  both  bending  and  membrane  capacity  was  chosen  because  it  provided  a  little 
additional  accuracy,  as  well  as  reduced  the  number  of  unrestrained  degrees  of  freedom 
which  had  to  be  constrained  to  avoid  singularities.  This  minimized  errors  caused  by 
manual  restraint  of  these  degrees  of  freedom.  In  some  applications,  such  as  lami¬ 
nated  plates  under  lateral  loads  where  only  bending  occurs,  the  use  of  CQUAD1  or 
CQDPLT  (with  bending  capacity  only)  elements  are  mandatory. 

The  NASTRAN  element  connection  card  CQUAD1  references  a  property  card  PQUADl 
which,  in  turn,  references  a  material  property  card  MAT2.  The  MAT2  cards  define  the 
material  properties  for  linear,  temperature-independent,  anisotropic  materials.  This 
card  requires  that  the  program  user  specify  the  elements  C„  of  the  two-dimensional 
anisotropic  stress-strain  relations  below  (see  Reference  11,  pp.  2.4-143): 
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The  elements  G..  of  the  stress-strain  relations  are  synonymous  with  the  variable 
_ £ 

Q  which  is  used  in  this  report  and  which  is  also  common  terminology  in  laminated 
ij  5 

plate  theory.  It  is  important  to  distinguish  the  reduced  stiffnesses  G,.  above,  as 
used  in  NASTRAN,  from  shear  moduli.  The  only  exception  is  that  the  reduced  stiff¬ 
ness  G33  (as  defined  in  NASTRAN,  but  not  in  this  report)  turns  out  to  be  the  shear 

modulus  G.„  for  a  unidirectional  lamina  (see  Appendix  A).  The  procedure  for  deter- 
•*■2  * 

mining  the  G,,  or  Q..  consists  of  the  following  steps: 

ij  ij  _ 

1.  The  elements  Q  of  the  two-dimensional  stress-strain  relations  (Equa¬ 

tions  (4)  given  earlier) 
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are  first  determined  for  the  individual  layers  of  the  laminate  using  the  procedure 
in  Appendix  A. 

2.  The  Q  from  step  (1)  are  next  used  in  Equations  (3)  to  calculate  the 
extensional  and  bending  stiffnesses  A^  and  for  the  full  laminate,  as  de¬ 
scribed  in  Appendix  B. 

*  * 

3.  Lastly,  reduced  stiffnesses  Q^.  ^  and  Q^.  ^  for  the  full  laminate  are 
found  from  the  A^  and  using  the  two  equations  (see  Appendix  E) 

Q*  ,  v  -  A, ,/h 
ij  (m)  ij 

and 


«uo»  ■  V(h3/12) 


(30) 


where  h  is  the  thickness  of  the  laminate. 


Following  the  procedure  outlined  above,  the  ^  and  Q^.  for  the  lami¬ 
nated  graphite  epoxy  skin  of  the  box  beam  were  found  to  be 


* 

^ij  (m) 


0.17002X108  0.44043xl07  0.18775xlo6 

0.44043X107  0.55951X107  0.18776xl06 

L0.18775X106  0.18770X106  0.47976xl0? 


and 


Qij  (b) 


0.16597X108  0.48865X107  -0.61273X101 
0.48865X107  0.60800X107  0.61273X101 

L-0.61273X101  0.61273X101  0.52632xl07 


(31) 


In  many  laminate  applications  the 


"k 

and  Q,. 


will  turn  out  to  be  similar  in 


lj(b) 

value,  as  was  the  case  for  the  box  beam  laminate.  Greater  differences  in  magnitude 
for  these  reduced  laminate  stiffnesses  may  be  expected  when  the  individual  layers  of 
a  laminate  have  more  widely  varying  stiffness  properties.  For  example,  if  a  laminate 
is  fabricated  of  rubber  and  steel  layers,  with  the  steel  placed  in  the  outer  layers 
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of  the  laminate,  the  Q../  .  and  Q, .  ....  will  exhibit  greater  differences.  Another 

ij  (m)  xij  (b)  % 

way  of  looking  at  it  is  that  the  D..  (and  therefore  Q,.,,.)  are  sensitive  to  the 

ij  ij  (b)  * 

stacking  sequence  of  the  laminas  while  the  A^.  (and  therefore  are  not. 


Other  Modeling  Considerations.  The  box  beam  skin  in  the  vicinity  of  its  fixed  end 
presents  special  modeling  considerations  because  of  a  scarf  joint  (see  detail  A  of 
Figure  1)  at  that  location,  which  involves  a  hybrid  cross  section  made  up  of  both 
composite  and  steel  materials.  This  was  solved  by  use  of  a  "weighting"  formula  that 
combines  proportionate  effects  from  both  steel  and  composite  for  a  cross  section  to 
arrive  at  the  membrane  and  bending  stress-strain  relationships.  At  the  opposite  free 
end  of  the  beam,  referred  to  as  the  load  block  end,  the  entire  load  block  region 
(steel  and  composite)  is  considered  equivalent  to  an  all  steel  block,  for  simplifi¬ 
cation  of  the  analysis,  since  the  major  part  is  steel.  This  load  block  was  modeled 
with  isoparametric  hexahedron  elements,  as  indicated  in  Figure  15. 

Finally,  some  comments  are  made  with  respect  to  boundary  conditions.  When  the 
box  beam  is  viewed  from  the  side  (see  Figure  1),  it  is  noticed  that  a  massive  metal 
cross  section  extends  a  few  inches  from  the  base  plate.  To  simplify  the  model,  fixed 
end  conditions  are  imposed  at  the  cross  section  where  this  steel  sleeve  section  ends 
and  the  spars  and  skin  continue.  This  eliminates  complexity  of  the  model  by  exclud¬ 
ing  a  section  that  would  experience  relatively  little  deformation  under  the  applied 
load.  Symmetry  of  the  beam  suggests  a  further  reduction  of  the  model  by  eliminating 
half  of  the  beam  split  lengthwise  along  its  longitudinal  axis.  This  artificial 
boundary  is  fixed  against  transverse  displacement  and  rotation  about  the  longitudinal 
axis  so  that  behavior  is  consistent  in  the  absence  of  the  other  half.  All  nodes  of 
the  structure  are  fixed  against  rotations  about  axes  normal  to  the  plate  elements 
because  these  elements  do  not  offer  resistance  to  this  type  of  displacement  and 
singularities  would  result  if  these  were  not  constrained.  The  isoparametric  hexa¬ 
hedron  elements  pose  a  similar  problem  to  the  plate  elements,  except  no  rotational 
degrees  of  freedom  are  allowed.  Therefore  the  nodes  bordering  these  elements  are 
fixed  for  rotation  about  all  of  the  space  axes. 


Technique  Effectiveness.  Using  the  NASTRAN  model  just  described,  a  finite-element 
analysis  was  performed  on  the  composite  box  beam  for  a  56-kip  end  load.  The  analysis 


results  were  useful  in  that  they  (a)  provided  a  check  on  the  strength-of-materlal 
calculation,  (b)  provided  detailed  results  on  the  structural  behavior  of  the  compos¬ 
ite  box  beam,  and  (c)  demonstrated  the  utility  of  the  NASTRAN  and  other  similar 
finite-element  computer  programs  for  the  analysis  of  hybrid  structures  involving 
composite  laminates. 

A  detailed  discussion  of  the  analysis  results  is  provided  elsewhere.*,**  The 
aim  of  this  report  is  to  present  the  finite-element  technique  for  laminated  com¬ 
posites  and  to  provide  a  demonstration  of  the  effectiveness  of  this  technique  for 
hybrid  structure  such  as  the  composite  box  beam. 

The  longitudinal  bending  stresses  in  the  composite  skin  and  steel  spar  of  the 
box  beam  are  presented  in  Figure  16  as  obtained  from  the  NASTRAN  analysis,  strength- 
of-materials  calculations,  and  static  tests.  Comparing  these  results  for  the  lami¬ 
nated  skin,  it  is  seen  that  the  NASTRAN  curve  is  in  very  good  agreement  with  the  test 
data.  Similar,  and  possibly  a  little  closer,  agreement  was  found  for  the  flange  of 
the  steel  spar.  The  strength-of-materials  curve  in  Figure  16  is  seen  to  underpredict 
the  test  data  by  approximately  5.5  percent  at  about  35.0  in.  from  the  fixed  end  where 
the  governing  stresses  were  obtained.  These  comparisons  demonstrate  the  effective¬ 
ness  of  the  finite-element  approach  described  earlier  for  composite  laminates  and 
also  show  that  conventional  strength-of-materials  calculations  are  effective  tools 
for  stress  prediction  in  laminates. 

Since  the  laminate  stresses  in  Figure  16  represent,  in  a  sense,  averaged  values 
or  equivalent  material  values  for  a  box  beam  skin  having  equivalent  anisotropic 
properties,  the  stresses  are  not  necessarily  the  maximum  values  which  might  occur 
in  any  of  the  layers.  For  example,  the  stresses  in  the  first  or  outer  layer  of  the 
skin  laminate  may  exceed  the  equivalent  material  or  averaged  stresses  for  the  entire 
laminate.  The  procedure  for  computing  the  stresses  in  any  layer  of  the  laminate  is 
given  in  Appendix  F. 

While  the  stresses  resulting  from  the  finite-element  analyses  were  only  averaged 
val|^  for  the  laminate  skin  and  required  further  processing  to  get  individual  layer 
sjtresses,  the  nodal  point  deflections  resulting  from  the  analysis  (see  Figure  17) 
represent  the  actual  deflections  and  need  no  further  processing. 

*Stein,  M.C.,  "A  Nastran  Analysis  of  a  Composite  Laminate  Box  Beam  for  Appli¬ 
cation  to  Navy  Hydrofoils,"  reported  informally  as  enclosure  (1)  to  DTNSRDC  ltr 
77-173-186  of  9  Dec  1977. 

**Barry,  M.  and  W.  Couch.  "Advanced  Composite  Box  Beam:  Static  and  NDE  Test 
and  Evaluation"  (in  preparation). 
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DISTANCE  FROM  BASE  PLATE  IIN.) 


Figure  16  -  Comparison  of  Longitudinal  Bending  Stresses  in 
Composite  Box  Beam  Based  on  Analyses  and  Tests 
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DISTANCE  FROM  BASE  PLATE  (IN.) 

Figure  17  -  Deflection  Curve  for  Composite  Box  Beam  with  a  56-Kip  Load 
Foil  Flap  Analysis 

The  finite-element  procedure  applied  to  the  composite  box  beam  and  discussed  in 
the  previous  section  was  also  used  to  analyze  the  composite  flap  for  the  aft  foil  of 
the  hydrofoil  PCH-1.  The  composite  flap,  illustrated  in  Figure  2,  consists  of  a 
laminated  skin  with  a  titanium  crank  spar  substructure  assembly.  The  composite  skin 
is  1/2  in.  thick  and  is  made  up  of  36  plies  of  a  graphite  epoxy  fabric  (T300) 

oriented  at  +45  deg  and  clad  by  10  mils  of  titanium. 

This  report  does  not  discuss  results  of  the  foil  flap  analyses  since  they  arc 
given  in  Reference  2.  Instead,  some  of  the  significant  differences  in  structural 

loading,  modeling,  and  response  of  the  laminated  skin  in  the  foil  flap  and  the  box 

beam  are  pointed  out  here.  Some  comments  are  first  made  with  respect  to  the  box 
beam.  Because  of  the  end  loading  on  the  box  beam,  the  laminates  in  the  upper  and 
lower  skins  or  flanges  are  subjected  to  longitudinal  strains  which  vary  only  slightl 


through  the  thickness  of  the  skin.  Since  the  in-plane  longitudinal  strains,  as  well 
as  the  transverse  and  shear  strains,  are  essentially  uniform  through  the  skin  thick¬ 
ness,  a  membrane- type  plate  element  CQDMEM  could  probably  have  been  used  in  the 
NASTRAN  analyses.  Instead,  the  CQUAD1  element  was  employed,  with  both  membrane  and 
bending  capacity,  to  improve  the  accuracy  slightly,  but  mainly  to  reduce  the  number 
of  degrees  of  freedom  which  had  to  be  otherwise  constrained,  as  mentioned  earlier. 

Although  the  CQUAD1  element  may  have  been  somewhat  optional  in  modeling  the 
box  beam  skin,  this  element  was  necessary  for  representing  the  foil  flap  skin  in 
order  to  accommodate  the  mostly  local  bending  response  associated  with  the  distrib¬ 
uted  pressure  loading  on  the  foil  flap. 

Another  difference  between  the  models  for  the  laminated  skin  of  the  foil  flap 


and  the  box  beam  has  to  do  with  the  comparative  magnitudes  of  the  reduced  stiffnesses 

Q  ,  „  and  0  „x.  The  greatest  difference  in  these  reduced  stiffnesses  for  the  box 

ij  (m)  xij  (b)  °  6  *  5 

beam  was  only  2.4  percent  =  17. 002*10  and  =  16.597*10  ).  With  so 

little  difference,  the  membrane  reduced  stiffness  could  have  been  input  for  both  the 

membrane  and  bending  components  of  the  CQUAD2  element  without  much  loss  in  accuracy . 

However,  in  the  case  of  the  laminated  foil  flap  skin,  the  greatest  difference  in 

magnitude  between  ^  and  Q„  ^  was  found  to  be  13.5  percent  =  7.506*10 


and  Q*  ,  ,  =  6.616*106).  Therefore,  it  is  more  important  for  the  foil  flap  analysis 
11 (m)  *  * 

than  it  was  for  the  box  beam  that  the  distinct  values  for  and  Qii be  used. 


<ij  (m) 


SUMMARY  AND  CONCLUSIONS 

Since  the  liklihood  exists  that  laminated  composites  will  be  increasingly  used 
in  the  fabrication  of  ship  and  marine  structures  in  general,  it  is  important  that 
appropriate  design  and  analytical  methods  be  available  for  supporting  their 
application. 

Relatively  recent  applications  of  laminated  composites  in  the  U.S.  Navy  have 
included  the  forward  foil  test  component  and  the  control  flap  on  the  foil  of  the 
hydrofoil  PCH-1.  It  is  possible  that  future  applications  could  eventually  include 
some  of  the  other  basic  components  of  ship  structure,  such  as  stiffeners,  girders, 
stanchions,  shell,  and  bulkhead  plating.  All  of  these  past  and  future  applications 
may  be  modeled  and  analyzed  as  one  of  the  laminated  structural  components  (i.e., 
laminated  beam,  column,  plate,  panels,  or  grillage.  .  .  see  Figure  4). 
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Three  types  of  analytical  methods  which  are  commonly  applied  to  metallic  struc¬ 
tures  have  been  extended  to  laminated  composite  applications.  These  include  (Table 
1):  strength-of -mater ials  techniques;  classical  methods  for  the  bending,  buckling, 
and  vibration  of  laminated  plates;  and  finite-element  and  finite-difference  methods. 
In  order  to  apply  these  methods  to  laminated  composites,  it  is  necessary  to  represent 
the  laminated  member  by  an  equivalent  nonlaminated  member  having  equivalent  elastic 
properties. 

Corresponding  to  each  of  the  three  types  of  analytical  methods  there  are  three 
kinds  of  elastic  properties  which  must  be  evaluated  in  advance  by  the  structural 
analyst.  The  required  stiffness  properties  are:  engineering  constants  E^,  E  , 
v  ,  V  ,  and  G  for  strength-of-materials  analyses;  extensional  and  bending  stiff- 

y  y  xy  * 

nesses  A.,  and  D..  for  classical  plate  analyses;  and  reduced  stiffnesses  Q..  for 
ij  ij  3  '  vij 

finite-element  analyses.  These  three  sets  of  stiffness  parameters  are  not  independ¬ 
ent  but  are  analytically  related  and  derived  from  each  other.  Extensional  and  bend¬ 
ing  stiffnesses  are  normally  calculated  first  and  then  used  to  compute  the  engineer¬ 
ing  constants  and  the  reduced  stiffnesses  as  needed.  Since  the  computation  of  these 
equivalent  stiffnesses  is  a  tedious  procedure,  if  performed  manually  for  a  laminate 
of  many  layers,  computer  programs  have  been  written  to  automate  the  procedure.  One 
nonproprietary  program  which  is  readily  available  to  the  public  is  SQ5  (see  Appendix 
D).  If  one  chooses  to  perform  these  calculations  manually  for  a  laminate  of  only  a 
few  layers,  or  if  one  wants  to  better  understand  the  theory  behind  the  stiffness 
calculations,  then  refer  to  Appendices  A,  B,  C,  and  E.  As  an  illustration  of  the 
manual  process  for  computing  engineering  constants,  calculations  are  given  in 
Appendix  C  for  the  laminated  skin  of  a  composite  box  beam  (for  application  to  the 
forward  foil  of  the  PCH-1  hydrofoil)  consisting  of  53  plies  of  T300  graphite  fibers 
and  32  plies  of  GY70  fibers.  These  constants  were  found  to  be  in  close  agreement 
with  experimentally  measured  engineering  constants,  as  indicated  in  Table  2.  This 
result  demonstrates  that  the  equivalent  engineering  constants  for  multilayered 
laminates  may  be  adequately  determined  analytically  without  resorting  to  more  costly 
experimental  means.  It  should  be  noted,  however,  that  it  is  common  practice  to  use 
experimental  measurements  to  determine  the  unidirectional  elastic  moduli  for  the 
constituent  layers  within  the  laminate  since  existing  analytical  procedures  such  as 
the  "law  of  mixtures"  are  not  always  of  sufficient  accuracy.  The  law  of  mixtures^ 
simply  computes  the  unidirectional  elastic  modulus  of  a  f iber- imbedded  layer  from 
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the  elastic  moduli  of  the  fiber  and  matrix  materials  and  the  percentage  composition 
of  each.  It  is  important  to  note  that  in  determining  the  equivalent  stiffnesses 
for  a  laminate,  a  distinction  is  made  between  equivalent  stiffness  properties  under 
extensional  and  bending  responses.  For  example,  the  engineering  constants  shown  in 
Table  2,  and  referred  to  above,  are  based  on  a  purely  extensional  response  of  the 
box  beam  laminate. 

As  mentioned  above,  three  types  of  analytical  methods  are  discussed  in  this 
report  for  analyzing  the  structural  response  of  hybrid  structures  involving  com¬ 
posite  laminates.  Solution  methods  are  discussed  for  analyzing  laminated  beams, 
columns,  and  rectangular  plates.  In  the  case  of  plates,  engineering  formulas  are 
readily  available  for  evaluating  the  deflection,  buckling,  and  vibration  behavior 
of  simply  supported  rectangular  plates  where  the  laminate  construction  consists  of 
a  cross-ply  arrangement  (Figure  7).  For  plates  having  angle-ply  layups,  the  analyst 
must  develop  desired  solutions  by  using  the  governing  differential  equations  and 
boundary  conditions. 

The  finite-element  method  is  the  most  general  tool  for  analyzing  the  structural 
response  of  laminates.  At  least  four  finite-element  approaches  are  available,  as 
illustrated  in  Figure  14,  for  investigating  laminates  depending  on  the  type  of  stress 
results  desired.  The  first  approach  involves  the  modeling  of  laminates  using  2D 
plate  elements  where  different  equivalent  stiffness  properties  are  allowed  for 
membrane  and  bending.  This  approach,  which  was  successfully  employed  and  validated 
on  the  composite  box  beam,  results  in  determination  of  in-plane  normal  end  shear 
stresses  within  the  laminate.  These  stress  results  correlated  well  with  test 
results  and  strength-of-materials  calculations.  The  finite-element  analysis  itself 
produces  equivalent  stresses  for  the  full  laminate,  but  these  may  be  further  pro¬ 
cessed  to  give  stresses  within  the  individual  layers  of  the  laminate.  If  a  knowl¬ 
edge  of  the  interlaminar  shear  and  through-the-thickness  normal  stresses  is  required, 
as  is  frequently  the  case  for  many  applications,  one  of  the  other  three  approaches 
in  Figure  14  must  be  implemented. 
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APPENDIX  A 

REDUCED  STIFFNESSES  FOR  INDIVIDUAL  LAMINAS 


A.l.  REDUCED  STIFFNESSES  FROM  UNIDIRECTIONAL  MODULI 

The  theory  below  applies  for  a  single  lamina  or  layer  of  a  composite  laminate 
having  fibers  imbedded  in  one  direction  only  (see  sketch). 

Once  the  unidirectional  elastic  constants  F.^,  E^,  V21’  an<*  ^12  are  ava^_ 

able  from  coupon  tests  for  a  given  individual  lamina,  the  two-dimensional  (2D) 
stress-strain  relations  may  be  obtained  in  terms  of  these  constants  from 


01 

i-H 

O' 

*12 

O 

_ 1 

(  1 

°2 

= 

*12 

Q22 

0 

2 

T12 

0 

0 

*66- 

(  Y12 

where  the  Q ,  the  reduced  stiffnesses,  are  given  by 


■n  i-vuvn 


U12  E2 


“21  E1 


12  1  -  vu  V21  1  -  v12  v21 


FIBERS 


and 


<22  1  -  V12  V21 


Q66  C12 


A. 2.  TRANSFORMATION  OF  REDUCED  STIFFNESSES  TO 
DIRECTION  OF  LOADING 

Since  the  individual  laminas  may  be  stacked  up  at  nonzero  (45  deg,  as  an 
example)  angles  to  the  direction  of  loading  to  achieve  desired  strength  and  stiffness 


i 
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for  the  overall  laminate,  it  is  necessary  that  the  2D  stress-relations  Equation  (A.l) 
for  the  laminas  having  fibers  oriented  at  some  angle  to  the  direction  of  loading  be 
transformed  to  the  direction  of  loading  on  the  laminate.  The  transformed  stress- 
strain  relations  have  the  form 


°x  r  Q11  Q12  Q16l  (Sc 

°y  =  ^12  ^22  ^26  |  Cy 

Txy  -  Q16  ^26  ^66  -  (  Yx 
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APPENDIX  B 

LAMINATE  EXTENSIONAL  AND  BENDING  STIFFNESSES 


As  was  discussed  earlier  in  the  main  text,  the  extensional  and  bending  stiff¬ 
nesses  A^j  and  D^. ,  respectively,  for  a  laminate  are  expressed  by  (assuming  a 
symmetric  layup  about  the  midplane  so  that  B^  =  0) 


f 

-h/2 


Qij  dz  = 


21  (Qij)k 


k=l 


and 


(B .  1) 


D.  . 


I 


h/2 


-h/2 


Qij  z2dz 


■  hZ  (5u 


v  ,  3  3 

k  Zk_Zk-l 


k=l 


Once  the  transformed  reduced  stiffnesses  Q„  have  been  found  using  the  proce¬ 
dure  of  Appendix  A,  the  task  of  finding  A^  and  D_  above  entails  carrying  out 
the  indicated  summations  involving  the  Q„  for  each  lamina  and  the  distances  to  the 
upper  and  lower  surfaces  of  the  laminas,  and  respectively  (see  sketch 

below) . 

The  task  of  computing  the  A.,  and  D..  using  Equations  (B.l)  for  a  laminate  with 
many  layers  is  a  tedious  one  and  well  suited  for  the  computer,  especially  the  D„  . 


DETERMINING  A. . 

ij 

The  summation  process  in  Equation  (B.l)  for  determining  A _  is  illustrated  below 
for  a  laminate  with  four  layers  having  symmetry  about  the  midplane. 


Aij  =  ^ijh  (zl'z0) 
+  (Q^)  2  (z2-zx) 

+  (Q±j ) 3  (z3~z2) 

+  (Qij) A  (z4-z3) 


(B.  2) 
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For  the  case  where  the  four  layers  are  of  equal  thickness  (h/4), 


zo  ■  -h/2 

z1  =  -h/4 

Z2  =  ° 

Zj  =  +h/4 

z4  =  +h/2 

Therefore, 


(B  .3) 


Aij  =  +  (Q^2h/A  +  <Qi;j>3h/4  +  (Q^y-M 


=  CCQ.j),  +  (Qy),  +  (Qij)3  +  (Q.j >4 ]  h/4 


(B.4) 


For  the  case  when  the  laminate  thicknesses  h^  through  h^  are  all  different, 


(B.5) 


However,  for  laminates  having  a  symmetric  layup  (such  that  the  coupling  stiffnesses 
B_  =  0,  as  assumed  earlier), 

h  «  h.  =  h  (Q,  .).  =  (Q..).  =  (Q..) 

1  4  a  ij  1  ij  4  ija 

(B.6) 

h2  =  h3  =  hb  (^ij*2  =  ^ij^  =  ^ij^b 

In  this  case. 


A..  =  2  (Q  .  . )  h  +  2  (Q  .  . ) ,  h. 
13  Hij  a  a  xij  b  b 


(B .  7  ) 
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DETERMINING  D 

The  bending  stiffness  is  given  by  the  following  summation  for  the  four- 

layer  laminate  previously  considered  where  the  four  layers  are  of  equal  thickness 
h/4 : 


3  (Qij>l  (^zi'zo)  +  3  ^1^2  (z2~zl) 


+  i  (QyJs  <z3-z2>  +T  <V*  (Z4-Z3> 


(B  .8) 


Substituting  z^  through  z^  in  terms  of  h  from  Equations  (B.3),  there  results 


»1,  ■  +  «i,>2  +  «H>3  +  7«iih>  155 


For  the  case  where  all  four  layers  are  of  unequal  thickness,  the  D  counterpart  of 
Equation  (B.5)  for  A„  may  be  written  down,  but  not  as  concisely  as  Equation  (B.5). 
Recall  from  before  that  the  absence  of  symmetry  in  the  layer  thicknesses  or  the 
Q.  about  the  midsurface  of  the  laminate  will  necessitate  the  computation  of  the 

ij 

B ^ .  from  Equation  (3)  since  they  are  no  longer  zero.  However,  in  most  structural 
designs,  the  laminates  are  designed  with  built-in  symmetry  to  reduce  the  complexity 
of  the  structural  analyses  needed  to  assess  structural  performance  and  integrity. 


APPENDIX  C 

NUMERICAL  EXAMPLE:  COMPUTATIONS  FOR  HYBRID  BOX  BEAM  LAMINATE 

DESCRIPTION  OF  HYBRID  LAMINATE 
•  Composition 

53  plies  of  T300  graphite  fibers  at  0°  to  loading 
direction  (ply  thickness  =  0.00516  in.) 

32  plies  of  GY70  graphite  fibers  at  +45°  to  loading 
direction  (ply  thickness  =  0.00692  in.) 

"  Stacking  Sequence 

0°,  (0°,  ±45°)g,  0°,  (+45°,  0°)8,  0° 

C.l.  EXTENSIONAL  STIFFNESSES 

Engineering  Constants  for  Unidirectionally  Reinforced  Lamina  (from  Reference  1) : 


Constant 

T300  Ply 

GY 70  Ply 

E1 

20.41X106  psi 

42.9X106  psi 

E2 

1.42xl06  psi 

0.89xl06  psi 

* 

V12 

0.181 

0.482 

V21 

0.0126 

0.010 

G12 

0.792X106  psi 

0.614xl06  psi 

^Computed  from  V-^  =  v2l  ^1^2' 
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Reduced  Stiffnesses  for  Unidirectionallv  Reinforced 

Lamina  (see  procedure 

Appendix  A,  Section  A.l): 

Stiffness 

T300  Ply 

GY70  Ply 

E1 

Q11  (l-v12  v21) 

20.457X106 

43.108X106 

Q  ,  VUC2 

Ql2  (1-V12  V21> 

0.258X106 

0.431X106 

E2 

,22  '  (1-h2  “21) 

1.423X106 

0.894xl06 

Q66  =  G12 

0.792X1Q6 

0.614xl06 

Transformation  of  Reduced  Stiffnesses  to  Direction  of  Loading  (see  procedure  in 
Appendix  A,  Section  A. 2): 


T300  Plies:  Since  T300  plies  are  already  oriented  in  the  direction  of  loading, 


no  transformation  of  the  Q„  for  T300  is  necessary. 


Therefore,  Q . .  =  Q. . . 

ij  ij 


GY70  Plies:  Substituting  the  Q„  values  from  above. 


and 


U1  =  (3Q11  +  3Q22  +  2Q12  +  4Q66)/8  =  16-91535xl°6 
U2  =  (Q11  "  ^22)/2  =  21-10674xl°6 

U3  =  (Q11  +  Q22  "  2Q12  "  4Q66)/8  =  5-08550xl°6 

U4  =  (Q11  +  Q22  +  6Q12  "  4Q66)/8  =  3-51656xl°6 

U5  =  (Q11  +  Q22  “  2Q12  +  4Q66)/8  =  5-69950xl°6 
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Then,  for  plies  with  fibers  at  0  *  +45  deg 


and 


ill  = 

Ul  + 

U2  cos20  + 

U3  cos40  =  11.83005*10 

12 

U4  ~ 

cos48  = 

10.60206*106 

'22  = 

U1  " 

02  cos20  + 

U3  cos40  =  11.83005*106 

16  = 

“IU2 

sin20  - 

sin40  =  -10. 55337*106 

'26  = 

“  1  U2 

sin20  + 

sin40  =  -10.55337*106 

'66  = 

U5  - 

cos40  = 

10.78500*10^ 

For  the  plies  having  fibers  at  0  =  -45  deg. 


and 


Qn  =  11.83005*10 


Q12  =  10.60206*10 


Q22  =  11.83005*10 


Q.,  =  10.55337*10 

1  O 


Q_,  =  10.55337*10 
26 


Q,,  =  10.78500*10 
bb 


Determination  of  Extensional  Stiffnesses  A..:  Extensional  stiffnesses  A.,  are  de- 

_ LL  iJ 

termined  using  Equation  (B.l)  (also  see  Equations  (B.2)  and  (B.3)  for  examples 

illustrating  the  use  of  Equation  (B.l))  as  follows: 


ij  ^1300(0°)  (T300)  QijGY70(+45°)  (+45  )  Ql3GY70(-45°)  (  45  } 
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Now, 


h(T3oo)  =  ^3  plies  x  0.00516  in. /ply  =  0.27348  in. 

and 

^(+45°)  =  h(-45°)  =  plies  x  0.00692  in. /ply  =  0.11072  in. 


Therefore, 


A..=Q  (0.27348)  +  Q. .  (0.11072)  +  Q. .  (0.11072) 

1J  JT300(0°)  1JGY70(+456)  1JGY70(-45°) 

and 

A  =  [(20. 457)(0. 27348)  +  (11.83005) (0. 11072)  +  (11 .83005) (0. 11072) ] *106 
=  8.21423X106 


Similarly, 

A12  =  2.41828x10^ 
A22  =  3.00881xl06 

A16  =  0,0 
a26  =  0.0 

and 

A,,  =  2.60483xl06 
bb 


C.2.  ENGINEERING  CONSTANTS  FOR  TOTAL  LAMINATE 

Inverse  of  Extensional  Stiffness  Matrix  A..:  Substituting  values  for  A.,  into 
_  _ ij_  ij 

Equation  (7)  gives 


A  =  A11(A22  A66  A26)  “  A12(A12  A66~A16  A26}  +  A16(A12  A26  A16  A22) 


=  49.1453x10 


18 
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Then, 


and 


A11  = 

(a22 

A66‘ 

■A26 

)/A 

=  1. 

59475x10" 

11 

CM 

\<* 

(A12 

A66- 

A16 

A26 

)/A 

=  -1.28175x10' 

A22 

(A11 

A66- 

A16 

>/A 

=  4. 

,35376xl0"7 

A16 

(A12 

A26' 

-A16 

A22 

)/A 

=  0 

A26  = 

(A11 

A26‘ 

A16 

A12 

)/ A 

=  0 

A66 

(All 

A22- 

"A12 

)/A 

=  3. 

83902xl0“7 

-7 


Determination  of  Engineering  Constants:  Inserting  the  values  of 
Equations  (9)  and  rounding  off  the  numerical  results  gives 


Ex  =  l/(Auh)  =  12.67x10  psi 


xy 


v 


xy 

yx 


=  l/(A22h)  = 

4.64X10 

psi 

6 

"  1/(A66h)  " 

5.26X10 

psi 

"A12/All  = 

0.804 

-a12/a22 

0.294 

n  =  A.  ,/ A. .  =  0 
xy,x  16  11 


n  =  A_ ,/ A_»  =  0 
xy , y  26  22 


H  =  A,  ,/A, ,  =  0 
x,xy  -6  66 


and 


n  =  A„,/A,,  =  0 
y ,xy  26  66 


63 


APPENDIX  D 

INPUT  DATA  FOR  COMPUTER  PROGRAM  SQ5 


The  general  content  of  each  card  (from  Reference  7)  in  the  SQ5  problem  deck 
is  as  follows: 


CARD 


COLUMNS 


FORMAT 


VARIABLE 


REMARKS 


1 

2-66 


Blank 

A65 


Title 


1-5 


6-10 


Key  1= 


Key  2= 


11-15  / 


15 


Key  3= 


Any  alphanumeric  information  de¬ 
scribing  the  problem  which  will 
be  printed  at  the  top  of  the 
first  page  of  the  problem  output. 

0 — Program  operation  continues 
after  computation  of  laminate 
data. 

1 — Program  terminates  after  com¬ 
puting  and  writing  out  the  ele¬ 
ments  of  the  constitutive 
matrices  and  the  average  lami¬ 
nate  properties. 

0 — No  point  stress  or  thermal 
analysis  will  be  done. 

1 — A  point  stress  analysis  will 
be  made  on  input  sets  of  Nx>  N^, 

and  N  .  One  card  per  load  case 
xy  r 

must  be  added  to  the  problem 

deck  (see  Card  6) .  This  key 

must  be  set  to  1  if  a  thermal 

analysis  is  to  be  performed. 

0 — No  additional  analysis  to  be 
performed . 

1 —  A  point  stress  analysis  will 
be  made  of  average  stresses  o  , 

0 0,  T  D,  and  9.  0  is  the  angle 

p  otp 

at  which  the  stresses  are  applied. 
This  analysis  is  for  in-plane 
loads  only. 

2 —  An  interaction  diagram  will  be 
computed  for  the  input  laminate. 
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CARD 


COLUMNS  FORMAT  VARIABLE 


16-20 


Key  4= 


21-25 


Key  5= 


26-30 

31-35 

36-40 


!MA 
NOMAT 
NCL 


1-5  15  KEYNU=  I 


REMARKS 


0 — No  thermal  analysis  will  be 
made . 

1 — Thermally  induced  in-plane 
T  T  T 

stress  N  ,  N  ,  N  and  moment 
x  y  xy 

M^,  M^,  M^  resultants  will  be 
x  y  xy 

computed  for  an  input  tempera¬ 
ture  change.  If  Key  4=1,  Key 
2  must  be  set  equal  to  1. 

0 — No  interlaminar  shear  stress 
analysis  will  be  made. 

1 — An  interlaminar  shear  stress 

analysis  will  be  made  for  input 

values  of  Q  and  Q  . 

x  y 

Number  of  lamina  (100  max.). 
Number  of  materials  (100  max.). 

Number  of  loading  cases.  This 

applies  to  sets  of  N  ,  N  ,  N 
rr  x  y  xy 

and  M  ,  M  ,  M  ,  temperature 

V  YU  1 


0 — Minor  Poisson's  ratio  of  each 
material  is  calculated  by  the 
program . 

1 — Minor  Poisson's  ratio  of  each 
material  is  input. 


1-9  \ 


10-18 


19-27 


El  (I) 


E2(I) 


U1(I) 


Modulus  of  elasticity  of  the 
I-th  material  along  the  first 
(or  "1")  lamina  axis. 

Modulus  of  elasticity  of  the  I-th 
material  along  the  second  (or  "2") 
lamina  axis  which  is  orthogonal  to 
the  "1"  lamina  axis. 

First  or  Major  Poisson's  ratio 
of  the  I-th  material. 
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CARD 

COLUMNS 

FORMAT 

VARIABLE 

REMARKS 

28-36 

G(I) 

Shear  modulus  of  elasticity  of  the 
I-th  material. 

37-45 

ALPHA  1(1) 

Coefficient  of  thermal  expansion 
of  the  I-th  material  in  the  "1" 
lamina  direction. 

4* 

46-54 

'  F9.0 

i 

ALPHA  2(1) 

Coefficient  of  thermal  expansion 
of  the  I-th  material  in  the  "2" 
lamina  direction. 

,  55-63 

ALPHA  6(1) 

Shearing  coefficient  of  thermal 
expansion  of  the  I-th  material. 

5** 

1-9 

F9.0 

U2(I) 

Second  or  minor  Poisson's  ratio 
of  the  I-th  material. 

1-5 

15 

LAY 

Lamina  number. 

6-10 

15 

MATYPE(I) 

Material  of  I-th  lamina. 

6***  < 

11-20 

F10. 0 

TH(I) 

Counterclockwise  angle  in  degrees 
from  the  laminate  reference  axes 
x,  y  to  the  lamina  natural  axes 

1,  2  of  the  I-th  lamina. 

21-30 

F10. 0 

AT(I) 

Thickness  of  the  I-th  lamina. 

1-10 

1  CALEl(I) 

Compression  limit  strain  allowable 
for  the  I-th  material  in  the  "1" 
lamina  direction. 

7  < 

11-20 

F10. 0  - 

CALE  2(1) 

Compression  limit  strain  allowable 
for  the  I-th  material  in  the  "2" 
lamina  direction. 

21-30  | 

CALE  3(1) 

Negative  limit  shear  strain  allow¬ 
able  for  the  I-th  material. 

31-40 

TALEl(I) 

Tension  limit  strain  allowable 
for  the  I-th  material  in  the  "1" 
lamina  direction. 

*There  will  be  a  Card  4  for  each  material. 

**There  will  be  a  Card  5  for  each  material  only  if  KEYNU  =  1  on  Card  3. 

***There  will  be  a  Card  6  for  each  lamina. 
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7* 


COLUMNS 

FORMAT 

VARIABLE 

REMARKS 

41-50 

TALE2(I) 

Tension  limit  strain  allowable  for 
the  I-th  material  in  the  "2" 
lamina  direction. 

51-60 

TALE3(I) 

Positive  limit  shear  strain  allow¬ 
able  for  the  I-th  material. 

1-9 

N(I,1) 

In-plane  force  resultant  in  the 
x  direction  for  load  case  I 
(Ib/in.) . 

10-18 

N(I,2) 

In-plane  force  resultant  in  the 
y  direction  for  load  case  I 
(lb/in.) . 

19-27 

N(I , 3) 

In-plane  shear  force  resultant 
for  load  case  I  (lb/in.). 

28-36 

F9.0 

M(I,1) 

M  moment  resultant  for  load 

X 

case  I  (in. -lb/ in.) . 

37-45 

M(I,2) 

M  moment  resultant  for  load 

y 

case  I  (in. -lb/in.) . 

46-54 

M(I,3) 

M  moment  resultant  for  load 
xy 

case  I  (in. -lb/in.) . 

55-63 

T(I) 

Change  in  temperature  for  load 
case  I. 

1-10 

SIG1 

Average  laminate  stress  acting 

in  a  direction  of  an  a,  3  system 
at  an  angle  PHI  from  the  laminate 
x,  y  axis  system. 

11-20 

Average  laminate  stress  crR  acting 
in  8  direction. 

*There  will  be  a  Card  7  for  each  material. 

**There  will  be  a  Card  8  for  each  load  case.  Card  8  is  omitted  if  (a)  only 
laminate  properties  are  desired,  (b)  only  an  interaction  diagram  is  desired,  or 
(c)  only  an  interlaminar  shear  analysis  is  desired. 
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CARD 

COLUMNS 

FORMAT 

VARIABLE 

REMARKS 

1 

9*  , 

21-30  j 
1 

F10. 0 

\  SIG3 

Average  laminate  shearing  stress 

0  o- 

osB 

31-40  , 

PHI 

Angle  in  degrees  from  the  a, 
system  to  the.  x,  y  system. 

8 

1-10 

QX(1) 

X  shear  force  resultant  for 
case  1  (lb/ in.) . 

load 

11-20 

QY(1) 

Y  shear  force  resultant  for 
case  1  (lb/in.) . 

load 

10** 

21-30 

►  F10.0 

QX(2) 

X  shear  force  resultant  for 
case  2  (lb/ in. ) . 

load 

31-40 

QY  (2) 

Y  shear  force  resultant  for 
case  2  (lb/ in. ) . 

load 

41-50 

QX(3) 

X  shear  force  resultant  for 
case  3  (lb/in.) . 

load 

51-60 

QY  (3) 

Y  shear  force  resultant  for 
case  3  (lb/in. ) . 

load 

*Card  9  is  input  only  if  Key  3=1  and  Key  1  =  Key  2  =  Key  4=0. 

**Enough  Card  10' s  must  be  included  to  cover  all  load  cases  if  Key  5=1. 
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APPENDIX  E 

REDUCED  STIFFNESSES  OF  LAMINATES  FOR  FINITE-ELEMENT  ANALYSES 


It  was  discussed  earlier  in  this  report  how  the  reduced  stiffnesses  0  for 

1J 

laminates  may  be  used  in  finite-element  programs  (such  as  NASTRAN)  to  represent  the 
two-dimensional  anisotropic  elastic  behavior  of  an  equivalent  nonlaminated  material. 

•ff 

This  appendix  outlines  the  steps  for  determining  these  Q..  (referred  to  as  G . .  on 
the  MAT2  card  of  NASTRAN). 


Consider  the  laminated  plate  element,  shown  on  the  left  in  the  sketch  below, 
having  unit  length  and  width  and  thickness  h.  The  laminate  may  have  any  number  of 
layers,  layer  thickness,  and  fiber  orientation  just  so  long  as  the  laminate  ends  up 
being  symmetric  about  the  midplane  (so  that  coupling  stiffnesses  B^.  between  exten¬ 
sion  and  bending  are  zero — see  Table  2  earlier  in  text).  The  laminated  plate  ele¬ 
ment  is  to  be  replaced  by  an  equivalent  (in  the  sense  that  it  has  the  same  exten- 
sional  and  bending  stiffnesses)  nonlaminated  plate. 


Laminated  Plate 


Equivalent  Nonlaminated  Plate 


For  the  symmetric  laminated  plate  element  above,  the  in-plane  forces  are  re¬ 
lated  to  in-plane  strains  by 


l  N 

1  x 

A11  A12  A16 

|  N 

1  y 

II 

A12  A22  A26 

1  N 

\  xy 

.A16  A26  A66_ 

(E.l) 
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and  the  moments  are  related  to  the  out-of-plane  curvatures  by 


"d  , 

D  „ 

3. 

i  k 

11 

12 

16 

I  x 

°12 

°22 

°26 

1 k 
)  y 

D  , 

D„, 

1  k 

16 

26 

66 

\  *y 

(E.2) 


Or,  using  condensed  matrix  notation, 


{N }  =  [ A ] { e }  and  {M}  =  [ D ] { k } 


(E.3) 


Likewise,  the  force-strain  and  moment-curvature  relations  on  the  laminated  plate  may 
be  written  as 


{N*}  =  [A*]{e*}  and  {M*l  =  [D*]{k*}  (E.4) 

*  * 

where  [A  ]  and  [D  ]  are  the  extenslonal  and  bending  stiffnesses  for  the  equivalent 
non lamina ted  plate.  Now,  if  for  the  same  applied  loads  on  the  laminated  and  equiv¬ 
alent  nonlaminated  plate,  that  is,  for 

(N)  =  {N*}  and  {M>  =  {M*} 

the  two  plates  have  the  same  strain  and  curvature 

{e}  =  {e*}  and  {k}  =  {k*} 

then  it  follows  that  the  plates  must  also  have  the  same  extensional  and  bending 
stiff  nesses 

[A]  =  [A*]  and  [D]  =  [D*]  (E.5) 

"k 

Next,  using  either  Equations  (3)  or  (B.l)  (both  are  the  same),  the  stiffnesses  [A  1 
* 

and  [D  ]  may  be  shown  to  be  given  by 


[A  ]  =  [Qjh 


[D*]  =  [Q*]h3/12 


(E.6) 


Therefore,  substituting  Equation  (E.6)  into  Equation  (E.5),  we  have 


‘  'A1  K 


(E.7) 


[Q*]  =  [D]  ~ ~ 

b  hJ/12 


Q*.  =  A . . / h 

<»)  1J 


(E.8) 


*  3 

Q.  .  =12  D. ,/h 

'ii  -l  -i 


•k  : k 

Using  Equation  (E.8),  the  reduced  stiffnesses  Q . .  and  Q..  may  be  straight¬ 


forwardly  computed  once  the  A„  and  are  found  either  by  hand  (see  Appendixes  A 

through  C)  for  laminates  having  a  few  layers  or  by  computer  using  the  program  SQ5 
for  laminates  of  many  layers  (see  Appendix  D) .  When  using  NASTRAN,  the  stiffnesses 
in  Equation  (E.8)  are  input  to  NASTRAN  as  the  variables  G„  on  the  MAT  2  cards,  as 
previously  indicated. 
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APPENDIX  F 

PROCEDURE  FOR  DETERMINING  TWO-DIMENSIONAL  IN-PLANE  LAMINA 
STRESSES  BY  POST-PROCESSING  FINITE-ELEMENT  OUTPUT 

The  following  procedure  is  used  to  determine  in-plane  lamina  stresses  by  post¬ 
processing  the  computer  finite-element  output. 

1.  Determine  the  equivalent  material  stress  on  the  upper  (u)  and  lower  (1) 
surface  of  the  element  from  computer  finite-element  analysis  output  (use  sign  con¬ 
vention  shown  in  sketch  on  right): 


2.  Decompose  the  surface  stresses  into  membrane  (m)  and  bending  (b)  component 


*Values  from  finite-element  output. 
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3.  Compute  equivalent  membrane  and  bending  strains: 


I 


* 

where  [C..]  =  [Q..]  (see  Appendix  E) . 

1J  b  b 


4.  Compute  total  strains  in  the  i-th  layer: 
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